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1.  Motivation  and  Objectives 

Reconfigurable  flight  control  systems  provide  significant  advantages  over  conventional  control  sys¬ 
tems,  in  particular  increased  probabilities  of  survival  after  failures  and  combat  damage.  Consider 
for  example  the  picture  of  a  fighter  aircraft  shown  in  Fig.  1.  The  shaded  areas  show  portions  of  the 
aircraft  that  might  be  missing  after  a  hypothetical  damage. 


Figure  1:  F-15  Aircraft  with  Hypothetical  Damage 

Such  an  incident  would  have  a  number  of  effects  on  the  dynamic  response  of  the  aircraft,  including: 
(1)  a  loss  of  the  balance  of  the  forces  and  moments  acting  on  the  aircraft,  resulting  in  a  divergent 
motion  if  the  remaining  control  surfaces  are  not  repositioned;  (2)  a  loss  of  the  stabilization  effect  of 
the  left  stabilator,  yielding  an  aircraft  with  lower  stability  margins;  (3)  a  reduction  of  the  ability  of 
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the  control  surfaces  to  generate  pitching  and  rolling  moments;  (4)  a  loss  of  symmetry  of  the  aircraft, 
so  that  a  symmetric  deflection  of  the  stabilators  produces  a  significant  rolling  moment,  in  addition 
to  the  pitching  moment.  If  the  control  law  is  fixed  and  optimized  for  the  unfailed  aircraft,  the 
damage  may  make  it  difficult  to  handle  the  aircraft,  if  not  impossible.  In  contrast,  a  reconfigurable 
control  system  may  be  designed  to  restore  the  stability  and  the  control  of  the  aircraft  by  learning 
about  its  dynamic  characteristics  in  real-time.  As  a  result,  steady  flight  may  be  maintained  and 
tracking  of  the  pilot  commands  restored. 

While  earlier  approaches  to  control  reconfiguration  required  the  failure  to  be  detected  and 
identified,  the  objective  of  this  project  was  to  develop  adaptive  algorithms  that  did  not  require 
failure  identification,  so  that  reconfiguration  could  be  achieved  for  a  large  number  of  failures  and 
damages  (including  some  which  may  not  have  been  predicted  at  the  design  stage).  A  specific 
emphasis  of  the  project  was  the  study  of  robustness  issues  related  to  unmodelled  dynamics  and 
actuator  limits.  The  results  of  the  project  have  been  reported  in  papers  listed  in  the  references. 
An  M.S.  thesis  was  also  completed,  which  is  included  in  Appendix.  This  report  is  an  executive 
summary  of  the  results  found  in  the  references. 

2.  Executive  Summary  of  Research  Accomplishments 
2.1  Multivariable  Adaptive  Control  Algorithms 

Multivariable  model  reference  adaptive  control  algorithms  were  developed  for  the  application  to 
reconfigurable  flight  control.  The  specific  algorithms  chosen  were  considered  attractive  because  of 
their  simplicity,  which  made  them  easy  to  implement  in  real-time,  and  useful  in  applications  with 
limited  computational  resources.  The  performance  of  the  algorithms  was  evaluated  in  simulations, 
using  a  highly  detailed  model  of  a  jet  aircraft,  similar  to  the  F-15  aircraft  shown  in  Fig.  1.  The 
model  included  a  nonlinear  aerodynamic  data  base,  models  of  actuator  dynamics  with  position  and 
rate  limits,  engine  dynamics,  aerodynamic  forces,  and  atmospheric  properties.  Trim  changes  and 
cross- couplings  between  the  longitudinal  and  lateral  motions,  both  significant  effects  of  failures, 
were  represented  in  the  model. 

The  ability  of  the  control  laws  to  maintain  steady  flight  and  to  restore  tracking  of  pilot  com¬ 
mands  after  unpredicted  failures  was  demonstrated  in  [1],  and  their  performance  over  a  wide  range 
of  flight  conditions  was  improved  during  the  course  of  this  project  [2].  Fig.  2  shows  some  typical 
simulation  results.  The  plots  show  the  responses  of  the  aircraft  to  a  pitch  rate  command.  The 
command  consists  in  steps  of  10  deg/s,  0  and  -10  deg/s.  The  damage  occurs  5  seconds  into  the 
simulation  and  consists  in  a  missing  left  stabilator  control  surface.  On  the  plots,  the  solid  lines 
are  the  responses  obtained  with  the  adaptive  control  law,  and  the  dashed  lines  are  the  responses 
obtained  without  reconfiguration,  that  is,  with  the  parameters  of  the  adaptive  control  law  frozen 
at  the  time  of  the  failure. 

The  plots  on  the  first  row  show  the  pitch  rate  response  (left)  and  the  roll  rate  response  (right) 
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Figure  2:  Responses  with  and  without  Adaptation  after  a  Missing  Stabilator 

(both  to  a  pitch  rate  command).  Without  reconfiguration,  the  desired  pitch  rate  is  not  achieved 
because  of  the  loss  of  pitching  moment,  and  the  roll  rate  is  very  large  because  of  the  coupling.  With 
reconfiguration,  the  required  pitch  rate  is  achieved,  and  the  coupling  in  roll  rate  is  much  reduced. 
The  plots  on  the  second  row  show  the  command  to  the  stabilators  (left)  and  the  command  to  the 
ailerons  (right).  The  stabilator  command  is  doubled  after  the  failure,  as  is  to  be  expected  since  only 
one  the  stabilators  is  working.  The  aileron  command,  which  was  zero  before  the  failure,  becomes 
much  larger  after  the  failure  in  order  to  compensate  for  the  rolling  moment  produced  by  the  single 
stabilator. 

The  plots  on  Fig.  3  show  the  estimates  of  two  significant  parameters  used  by  the  adaptive  control 
law.  The  parameter  on  the  left  is  the  amount  of  pitching  acceleration  generated  by  a  unit  deflection 
of  the  stabilators.  The  parameter  is  negative  and  decreases  in  magnitude  by  a  factor  of  2  after  the 
damage,  as  expected.  The  parameter  on  the  right  is  the  amount  of  pitching  acceleration  generated 
by  a  unit  change  in  angle  of  attack  and  is  indicative  of  the  inherent  stability  of  the  aircraft.  A  large 
negative  value  indicates  a  strong  tendency  of  the  aircraft  to  return  to  level  flight  after  a  disturbance. 
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Figure  3;  Adaptive  Parameter  Responses  after  a  Missing  Stabilator 

The  magnitude  of  the  estimate  is  reduced  after  the  damage,  as  expected  because  of  the  loss  of  the 
stabilization  effect  of  the  stabilator.  Note  that  the  parameters  converge  within  15  seconds,  although 
the  values  that  are  reached  within  a  much  shorter  time  are  sufficient  to  considerably  improve  the 
responses.  Overall,  the  responses  of  Figs.  2  and  3  show  that  the  adaptive  system  is  able  to  correctly 
identify  and  compensate  for  changes  in  system  dynamics  caused  by  the  damage. 

2.2  Parameter  Identification 

At  the  core  of  the  reconfigurable  control  law  discussed  in  section  2.1  was  a  parameter  identification 
algorithm  capable  of  determining  the  necessary  aircraft  parameters.  The  problem  of  parameter 
identification  for  reconfigurable  flight  control  is  similar  to  the  off-line  identification  problem,  but 
involves  additional  difficulties  which  may  be  categorized  as  follows: 

•  Lack  of  control  over  the  actuator  signals.  A  significant  difficulty  arises  because  the  actuator 
signals  are  determined  by  the  control  law  and  cannot  be  freely  selected.  Considerable  im¬ 
provements  in  identification  performance  can  be  obtained  through  optimization  of  the  signals 
applied  to  the  control  surfaces.  Conversely,  in  a  reconfigurable  control  application,  perfor¬ 
mance  may  be  poor  because  the  signals  exhibit  highly  undesirable  characteristics,  including: 

1.  high  levels  of  correlation  between  the  control  signals  and  the  aircraft  states  (in  particular 
for  linear  state  feedback  control  laws); 

2.  loiig  periods  of  quiescence  {e.g.,  in  steady  flight); 

3.  highly  coupled  longitudinal  and  lateral  motions  (in  particular  for  rapid  roll  maneuvers); 

4.  excitation  of  nonlinear  dynamics  (for  example  inertial  couplings). 
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In  addition,  because  the  control  signals  are  determined  by  a  control  law  which,  itself,  de¬ 
pends  on  the  identification  results,  adverse  interactions  are  sometimes  induced  between  the 
identification  and  the  control  components  of  the  system,  further  aggravating  the  problem. 

•  Increased  number  of  parameters  for  failed  aircrafts.  Aircrafts  often  loose  their  symmetry  after 
failures,  requiring  that  more  parameters  be  identified.  For  example,  the  roll  effectiveness  of  a 
symmetric  tail  deflection  may  be  assumed  to  be  zero  for  an  unfailed  aircraft,  but  not  otherwise. 
In  general,  accounting  for  possible  failures  may  prevent  the  ganging  of  some  surfaces,  such  as 
spoilers,  and  result  in  a  significant  increase  in  the  number  of  parameters  to  be  identified. 

•  Real-time  operation.  The  computational  requirements  of  a  parameter  identification  scheme  for 
reconfigurable  control  must  be  compatible  with  available  computers.  Further,  the  procedure 
must  be  performed  with  minimal  supervision  and  may  not  require  extensive  trial  and  error 
adjustments  by  the  pilot. 

An  interesting  adaptive  algorithm  was  proposed  by  Ward  &  Barron  ^  to  identify  aircraft  param¬ 
eters  in  the  face  of  these  challenges.  Their  so-called  modified  least-squares  algorithm  incorporated 
a  penalty  on  parameter  variations  in  the  optimization  criterion.  Ward  &  Barron  observed  that, 
without  the  penalty,  an  unacceptable  trade-off  was  encountered  in  the  design  of  a  reconfigurable 
flight  control  system:  either  the  responses  of  the  adaptive  algorithm  were  too  slow  during  periods 
of  significant  excitation,  or  too  noisy  during  periods  of  low  excitation. 

A  recursive  foriri  of  the  algorithm  was  developed  and  analyzed  by  the  principal  investigator 
[3].  An  averaging  analysis  showed  that  the  algorithm  exhibited  some  very  interesting  properties, 
including:  (a)  a  variable  data  forgetting  capability,  so  that  the  memory  of  the  algorithm  was  long 
when  the  excitation  was  poor,  and  short  when  the  excitation  was  significant;  (b)  a  second-order 
filtering  capability,  as  opposed  to  the  first-order  filtering  of  the  original  algorithm,  reducing  the 
rate  of  fluctuations  of  the  parameters;  (c)  a  stabilizing  effect  on  the  covariance  matrix  update. 

Fig.  4  shows  the  results  of  a  simulation  with  a  single  parameter  system,  illustrating  the  two 
features  mentioned  above.  On  the  left  is  the  response  of  the  modified  algorithm.  On  the  right  is  the 
response  of  the  standard  least-squares  algorithm  with  forgetting  factor.  The  true  parameter  varies 
between  1  and  —1  and  appears  as  a  triple  step  function.  The  smooth  curves  are  the  responses 
predicted  by  the  averaging  analysis,  while  the  noisy  curves  are  the  responses  of  the  parameter 
estimate.  The  approximate  responses  provided  by  the  averaging  analysis  are  found  to  be  close 
to  the  original  responses,  validating  the  analysis.  At  time  400,  there  is  a  drop  in  excitation  (and 
therefore  of  the  signal-to-noise  ratio).  The  response  of  the  modified  algorithm  is  seen  to  slow  down 
after  that,  while  the  parameter  fluctuations  remain  reasonable.  For  the  standard  least-squares 
with  (the  same)  forgetting  factor,  it  is  found  that  when  the  excitation  drops,  the  fluctuations 
in  the  parameters  become  very  large.  These  fluctuations  can  be  reduced  by  using  a  different 

^D.G.  Ward  Sz  R.L.  Barron,  Self- Designing  Flight  Control  Using  Modified  Sequential  Least-Squares  Parameter 
Estimation  and  Optimal  Receding  Horizon  Control  Laws,  Technical  Report,  Barron  Assoc.,  Charlottesville,  VA, 
March  1994. 
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Figure  4:  Adaptive  Parameter  Responses  -  Left:  Modified  LS;  Right:  Original  LS 

forgetting  factor,  but  the  response  in  the  first  part  of  the  simulation  becomes  slower  then.  With 
the  modified  algorithm,  the  algorithm  reacts  fast  when  there  is  sufficient  excitation,  yet  keeps 
parameter  fluctuations  reasonable  when  the  information  available  cannot  sustain  a  fast  adaptation. 

The  identification  algorithm  was  used  in  the  simulations  discussed  in  section  2.1,  as  well  as  in 
experiments.  During  the  summer  of  1996,  a  series  of  flight  tests  were  performed  with  the  VISTA 
F-IG  aircraft,  using  the  modified  least-squares  algorithm  and  a  reconfigurable  control  law.  These 
experiments  were  carried  out  by  Barron  AssociateSj  Inc,  under  contract  to  the  Air  Force,  and 
were  not  part  of  this  project.  However,  the  principal  investigator  assisted  in  the  planning  of  the 
experiments  and  in  the  analysis  of  the  data.  The  flight  tests  culminated  in  a  successful  landing  of 
the  aircraft  under  a  simulated  missing  elevon,  with  an  adaptive  system  reconfiguring  rapidly  after 
single  or  multiple  impairments.  The  design  of  the  parameter  identification  procedure  as  well  as  the 
experimental  results,  were  reported  in  [4]. 

The  adaptive  algorithm  was  also  tested  experimentally  under  more  modest  circumstances  in  the 
principal  investigator  laboratory.  In  this  case,  the  control  law  was  a  control  law  designed  to  avoid 
the  excitation  of  resonant  modes  in  a  flexible  beam.  Such  problem  is  more  relevant  to  the  control  of 
flexible  modes  in  aircraft  structures  than  to  the  flight  control  problem.  However,  the  experiments 
demonstrated  again  the  effectiveness  of  the  identification  algorithm  in  a  practical  environment.  The 
results  were  reported  in  [5],  [6]. 

2.3  Integral  Compensation  in  Adaptive  Control 

The  incorporation  of  integral  compensation  in  the  reconfigurable  control  laws  was  considered,  with 
the  objective  of  improving  performance,  and  of  determining  whether  certain  problems  might  be 
encountered.  Indeed,  it  Wcis  found  that  the  combined  use  of  integral  compensation  and  adaptive 
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Figure  5:  Integral  Compensation 

bias  cancellation  could  sometimes  result  in  algorithms  with  uncontrollable  modes  [7].  Surprisingly, 
different  conclusions  were  reached  depending  on  the  type  of  algorithm  used.  Specifically,  the  algo¬ 
rithms  studied  in  this  project  (as  well  as  those  implemented  in  the  Vista  F-16  flight  tests)  were  found 
to  belong  to  a  class  of  algorithms  for  which  integral  compensation  and  adaptive  bias  cancellation 
could  be  combined  without  uncontrollability  problems.  Given  that  integral  compensation  could  be 
included  in  the  model  reference  control  algorithms,  a  question  was  then  whether  significant  benefits 
could  be  derived  from  it.  From  extensive  simulations,  it  was  found  that  there  was  an  advantage 
to  the  use  of  integral  compensation,  but  that  the  improvement  was  less  than  might  be  expected 
from  conventional  control  systems.  To  illustrate  this  point.  Fig.  5  shows  the  results  of  simulations 
for  control  laws  with  and  without  integral  action.  The  reference  input  was  a  sequence  of  steps  in 
pitch  rate  (first  row)  and  roll  rate  (second  row).  On  the  left  are  the  pitch  rate  responses,  and  on 
the  right  are  the  roll  rate  responses.  The  solid  lines  are  the  desired  responses,  the  dashed  lines  are 
the  responses  obtained  without  integral  action,  and  the  dot-dashed  lines  are  those  obtained  with 
integral  action.  A  failure  occurred  at  t=310s,  which  consisted  in  a  locked  stabilator. 

Overall,  tracking  of  the  commands  is  seen  to  be  better  when  integral  compensation  is  used. 
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However,  tracking  is  good  even  without  the  integral  compensation.  Trimming  of  the  aircraft  is 
also  achieved  without  integral  action  (there  is  no  visible  drift  in  the  responses).  An  explanation, 
provided  in  [7],  is  that  the  adaptive  algorithm  is  able  to  fine-tune  the  estimates  of  the  gain  matrices 
and  of  the  bias  term  in  order  to  provide  a  compensation  equivalent  to  that  of  integral  control. 
Overall,  it  was  concluded  that,  with  a  properly  designed  adaptive  algorithm,  integral  compensation 
could  be  incorporated  without  uncontrollability  problems,  but  also  without  benefits  that  were 
necessarily  significant  enough  to  warrant  its  use.  While  integral  compensation  is  often  included 
in  the  pitch  axis  of  non-adaptive  flight  control  systems,  its  purpose  is  mainly  to  assist  with  the 
trimming  of  the  aircraft,  so  that  an  adaptive  control  law  with  trim  bias  estimation  may  not  need 
integral  compensation. 

2A  Actuator  Saturation 

Actuator  limits  pose  a  major  problem  in  flight  control  system  design.  The  positions  of  control 
surfaces  are  limited  both  in  their  range  and  in  their  rate  of  motion.  If  not  accounted  for  in 
the  design,  these  constraints  may  lead  in  flight  to  a  significant  degradation  of  performance  and, 
sometimes,  to  pilot-induced-oscillations.  The  origin  of  the  problems  is  in  the  delay  of  the  response 
of  the  aircraft  resulting  from  the  actuator  saturation  as  well  as  in  the  multivariable  nature  of  the 
control  problem,  which  is  such  that  saturation  in  one  axis  may  lead  to  poor  responses  in  other  axes. 

For  reconfigurable  flight  control  systems,  actuator  saturation  becomes  even  more  problematic. 
Not  only  is  control  power  reduced,  but  the  couplings  between  longitudinal  and  lateral  axes  also 
become  significantly  larger  after  most  failures.  Simulations  of  the  early  model  reference  control 
algorithms  yielded  results  that  deteriorated  significantly  for  large  commands.  The  first  row  of 
Fig.  6  shows  responses  exhibiting  this  problem.  On  the  left  is  the  pitch  rate  response  to  a  pitch 
rate  command,  and  on  the  right  is  the  roll  rate  response  to  a  pitch  rate  command.  The  solid  line 
is  the  desired  response  and  the  dashed  line  is  the  actual  response.  In  the  simulations,  the  left 
horizontal  tail  surface  is  locked  at  time  t  =  ZIO  sec.  The  responses  show  a  very  large  roll  rate 
response  to  the  pitch  rate  command.  The  large  coupling  is  not  observed  for  smaller  commands  of  1 
deg/s  and  its  appearance  is  due  to  the  fact  that  the  various  actuators  saturate  independently,  and 
in  a  way  that  prevents  decoupling  using  a  linear  controller.  The  pitch  rate  response  also  exhibits 
glitches,  which  may  be  related  to  the  sharp  roll  rate  transients. 

Four  methods  for  command  limiting  were  studied  in  [8]  for  their  application  to  reconfigurable 
flight  control.  Significantly  different  concepts  were  considered:  the  scaling  of  the  control  inputs, 
the  relaxation  of  the  control  requirements,  the  scaling  of  the  reference  inputs,  and  the  least-squares 
approximation  of  the  commanded  accelerations.  On  the  second  row  of  Fig.  6,  responses  are  shown 
which  were  obtained  with  the  first  two  methods  of  command  limiting.  Plots  for  the  other  two  meth¬ 
ods  were  not  very  different  and  are  not  shown.  One  sees  that  tracking  and  decoupling  performance 
is  much  improved  with  the  methods  of  command  limiting.  All  the  methods  gave  comparable  results, 
which  was  surprising  since  some  methods  were  much  simpler  than  others,  and  some  methods  used 
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Figure  6:  Responses  without  command  limiting  (top)  and  with  command  limiting  (bottom) 

the  knowledge  of  the  generated  moments,  while  others  did  not.  More  research  will  be  needed  to 
confirm  whether  relatively  simple  methods  of  command  limiting  may  indeed  be  adequate  for  most 
purposes,  or  whether  the  more  sophisticated  methods  are  truly  needed. 

A  different  set  of  results  was  obtained  regarding  the  control  of  systems  with  input  and  state 
constraints  [9].  The  problem  under  consideration  was  the  determination  of  the  sets  of  states  of  a 
linear  system  with  control  limits  for  which  recovery  to  the  equilibrium  is  possible.  These  results  have 
not  been  applied  to  reconfigurable  control  problems  yet.  However,  it  is  possible  that  tailless  aircrafts 
with  unstable  modes  in  two  axes  and  with  limited  control  authority  will  have  to  be  controlled  with 
concerns  of  recoverability  in  mind.  A  procedure  was  developed  to  rapidly  calculate  recoverable  sets 
numerically.  The  method  was  limited  to  two-dimensional,  continuous-time  systems  with  a  single 
input,  but  it  was  found  to  provide  significant  computational  savings  over  other  methods. 
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Abstract 


The  problem  of  reconfigurable  flight  control  is  investigated,  focusing  on  model  reference  adaptive 
control  and  a  stabilized  recursive  least-squares  algorithm.  Previous  work  has  demonstrated  this 
controller  to  be  capable  of  maintaining  control  of  an  aircraft  despite  actuator  failures.  However, 
extensive  simulations  have  also  shown  the  need  for  further  developments  of  the  control  law  in 
order  to  improve  transient  responses  and  robustness  to  unmodeled  dynamics.  Modifications  to 
the  algorithm  are  proposed  to  achieve  those  objectives.  The  effectiveness  of  the  modifications  is 
evaluated  using  a  detailed  jet  aircraft  simulation  model. 


Chapter  1 

Introduction 


A  need  for  increased  reliability  and  for  cost  reduction  has  encouraged  the  development  of  reconfig- 
urable  flight  control  laws.  The  goal  of  reconfigurable  flight  control  is  to  adjust  the  commands  to  the 
aircraft  surfaces  automatically  following  structural  damage  or  failure  in  order  to  maintain  safety 
and  performance.  Increased  reliability  is  achieved  by  exploiting  the  computational  capabilities  of 
modern  computers  and  the  redundancy  of  multiple  control  surfaces.  Costs  are  also  reduced  by 
alleviating  the  need  for  back-up  systems. 

Reconfigurable  control  can  be  applied  in  numerous  failure  situations.  An  example  is  when  one 
of  the  stabilators  of  an  airplane  becomes  jammed.  The  controls  to  the  other  unimpeded  surfaces, 
including  the  other  stabilator,  the  ailerons,  and  the  rudder,  are  modified  to  compensate  for  the  loss 
of  control  power  in  the  failed  actuator.  Maneuverability  of  the  aircraft  is  restored  to  acceptable 
levels. 

The  approach  adopted  in  this  thesis  is  based  on  adaptive  control  theory.  Adaptive  control 
provides  two  major  benefits.  First,  the  failure  does  not  have  to  be  identified  among  a  set  of 
predetermined  possibilities.  This  feature  is  desirable  due  to  the  numerous  failures  that  can  occur. 
Second,  adaptive  controllers  can  continuously  adapt  to  varying  flight  conditions.  Some  related 
approaches  that  have  been  proposed  so  far  are  now  reviewed. 


1.1  Brief  Literature  Review 

1.1.1  Proportional-Integral  Implicit  Model  Following 

Huang  and  Stengel  [1]  implement  a  controller  using  proportional-integral  implicit  model  following. 
The  controller  is  based  on  a  linear  quadratic  (LQ)  control  approach  which  utilizes  a  cost  function 
that  penalizes  the  states,  the  inputs,  and  the  couplings  between  the  states  and  the  input.  Based  on 
the  output  error  integral,  an  LQ  controller  equation  is  developed  which  provides  both  a  proportional 
and  integral  feedback  of  error. 

Huang  and  Stengel  state  the  benefits  of  linear  quadratic  control  theory  as  (1)  a  direct  treatment 
of  multiple-input,  multiple-output  systems,  (2)  an  easy  transfer  from  continuous  to  discrete-time 
descriptions,  and  (3)  guaranteed  stability  when  certain  criteria  are  met.  The  application  of  im¬ 
plicit  model  following  allows  a  relatively  easy  implementation  with  low  gains  and  a  low  number  of 
parameters.  The  simulations  in  the  paper  considered  multiple  failure  types,  but  did  not  consider 
any  cross-coupling  effects  between  longitudinal  and  lateral  dynamics.  When  an  aircraft  sustains 
damage,  symmetry  is  usually  lost,  and  undesirable  couplings  appear.  A  disadvantage  of  the  linear 
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quadratic  method  is  that  it  requires  penalty  matrices.  The  assignment  of  penalties  is  usually  ad-hoc 
and  may  need  to  be  altered  as  a  result  of  a  failure. 

1.1.2  Linear  Regression 

Chandler,  Pachter,  and  Mears  [2]  utilize  linear  regression  as  an  identification  technique.  The 
algorithm  is  a  batch  type  algorithm  able  to  handle  time- varying  parameters  and  parameters  subject 
to  change.  A  minimum  variance  estimate  is  derived  using  measured  states  of  the  aircraft.  The 
authors  state  their  objectives  as  (1)  develop  an  efficient,  static  system  identification  method  for 
adaptive/reconfigurable  control,  and  (2)  significantly  improve  parameter  estimates  by  including  a 
priori  parameter  relationship  information  from  flight  mechanics. 

The  paper  does  not  specifically  demonstrate  control  performance.  However,  the  identification 
is  shown  to  correctly  determine  the  parameters  after  a  failure  in  the  aircraft.  The  identification 
algorithm  is  designed  to  operate  continuously  and  to  make  use  of  prior  information  to  simplify  the 
time  and  amount  of  computation  required.  A  linearized  model  is  used  for  the  flight  dynamics,  and 
cross-couplings  are  neglected.  The  identification  technique  is  shown  to  work  well  and  is  similar  to 
the  one  used  in  this  research. 

1.1.3  Feedback  Linearization 

Feedback  linearization  with  combined  feedback  and  feedforward  control  for  large  transport  aircraft 
is  investigated  by  Ochi  and  Kanai  [3],  [4].  Both  feedback  and  feedforward  control  are  used  because 
of  the  varying  speed  of  the  control  effectors:  the  feedforward  control  is  applied  to  such  compo¬ 
nents  as  the  flaps  and  engines;  feedback  control  is  applied  to  the  elevators,  ailerons,  and  rudder. 
Feedback  linearization  uses  a  standard  state-space  model  with  the  state-space  matrices  being  linear 
combinations  of  constant  parameters  and  known  functions  of  the  inputs  to  the  system.  The  output 
matrix  is  comprised  of  nonlinear  functions.  The  actual  output  must  track  the  desired  output.  The 
problem  is  simplified  by  utilizing  generic  inputs  which  aJlow  combining  symmetric  surfaces  under 
the  same  command.  The  feedback  linearization  technique  is  applied  to  a  large  transport  aircraft 
as  opposed  to  a  fighter  aircraft.  A  nonlinear  flight  model  is  used,  and  the  controller  is  found  to  be 
able  to  maintain  flight.  Cross-couplings,  however,  are  not  considered. 

1.1.4  Neurocontrollers 

A  somewhat  different  approach  to  aircraft  control  is  the  use  of  neurocontrollers.  This  technique  is 
proposed  by  Napolitano  et  al.  [5]  The  fundamental  concept  of  neurocontrollers  is  that  they  “learn” 
the  flight  conditions  of  the  unfailed  aircraft  and  apply  the  “learned”  information  to  bring  the 
aircraft  back  to  equilibrium  once  the  failure  occurs.  Napolitano  et  al.  state  this  technique  adapts 
well  to  nonlinear  behavior.  The  learning  is  based  on  the  extended  back  propagation  algorithm 
(EBPA)  which  avoids  the  pitfalls  of  the  standard  back  propagation  algorithm  (SBPA).  Failure 
identification  is  necessary  with  this  controller.  Once  a  failure  is  detected,  the  controller  switches 
from  the  “learning”  mode  to  the  “controlling”  mode.  Control  after  failure  is  based  upon  following 
a  model.  A  cost  function  is  implemented  in  this  controller.  In  the  unfailed  case,  cost  is  based 
on  following  the  actual  surface  control  deflections.  In  the  failed  case,  cost  is  based  on  the  errors 
between  the  actual  and  desired  axis  angular  position,  rate,  and  acceleration. 

The  use  of  neurocontrollers  for  flight  control  is  a  recent  development.  Control  of  the  nonlinear 
system  is  demonstrated  with  a  feasible  on-line  implementation.  Cross-coupling  effects  are  consid- 
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ered  with  this  controller.  One  issue  that  should  raise  some  concern,  however,  is  the  large  transients 
that  appear  in  the  simulations  after  a  failure  occurs.  Another  issue  that  should  raise  some  concern 
is  that  of  “learning  time.”  Napolitano  et  al.  state  the  EBPA  algorithm  has  a  shorter  learning  time 
than  the  SBPA  algorithm.  However,  any  learning  time  could  cause  the  aircraft  to  become  unre¬ 
coverable  if  a  failure  occurs  before  the  algorithm  acquires  sufficient  information.  The  simulations 
in  the  paper  showed  a  “learning  time”  of  about  20  seconds  which  is  quite  long. 

1.2  Approach  of  the  Thesis 

Model  reference  adaptive  control  is  considered  for  solving  the  reconfigurable  flight  control  problem. 
Benefits  of  the  algorithm  include 

•  the  simplicity  of  the  algorithm, 

•  the  capability  to  handle  cross-couplings,  and 

•  the  state  feedback  nature  of  the  control  law. 

Flight  control  algorithms  should  be  simple.  The  computational  capability  of  on-board  computers 
is  limited,  and  calculations  should  be  kept  to  a  minimum.  In  addition,  keeping  an  algorithm  simple 
reduces  the  amount  of  difficulty  in  trouble-shooting.  Cross- couplings  also  are  significant  in  flight 
control  for  failed  aircrafts.  An  unfailed  aircraft  is  normally  symmetric,  resulting  in  a  decoupling 
between  longitudinal  and  lateral  variables.  However,  in  times  of  failure,  the  symmetry  is  usually 
lost.  The  effects  of  cross-couplings  require  a  multivariable  control  law.  Finally,  the  state-feedback 
nature  of  the  control  law  is  well-suited  for  the  flight  control  problem  because  of  the  availability 
of  state  measurements.  Stability  is  also  better  understood  for  state-feedback  control  than,  for 
example,  for  a  predictive  control  law. 


1.3  Review  of  Previous  Work 

Work  previous  to  this  thesis  was  performed  in  [6].  There  were  three  main  goals  for  the  system: 

•  selection  of  proper  trim  values  for  the  inputs,  even  if  a  failure  caused  rapid  changes  in  those 
values, 

•  decoupling  of  the  inputs  and  outputs,  and 

•  tracking  of  the  reference  inputs. 

The  realization  of  these  three  goals  was  achieved  in  [6].  However,  the  set  of  conditions  for  which 
the  performance  of  the  algorithm  was  evaluated  was  very  small.  Therefore,  an  original  objective 
of  this  thesis  was  to  expand  the  evaluation  of  the  adaptive  algorithms  for  a  broader  range  of  flight 
conditions  and  to  more  demanding  maneuvers.  In  the  process,  unstable  behavior  of  the  closed-loop 
system  was  discovered  in  several  instances.  The  problem  was  complicated  by  the  fact  that  no 
single  cause  could  be  identified  for  the  instability.  The  following  factors  were  determined  to  be 
contributing: 

•  low  excitation:  the  parameters  of  the  control  law  were  found  to  diverge  during  times  when 
excitation  was  not  provided  in  certain  axes.  The  parameters  were  not  accurately  identified, 
resulting  in  inadequate  control  of  the  aircraft.  If  excitation  was  provided  to  all  the  inputs, 
the  instabilities  disappeared. 
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•  actuator-  dynamics:  the  adaptive  control  law  assumed  that  the  control  surfaces  moved  instan¬ 
taneously.  However,  the  actuators  of  the  aircraft  did  not  respond  instantaneously,  and  the 
time  delay  of  the  actuators  was  found  to  contribute  to  the  instabilities. 

•  actuator  limits  and  rate  limits:  the  control  law  calculated  control  inputs  without  constraints 
on  the  input  values.  However,  with  the  actuators  limited  in  the  rate  and  range  of  motion, 
unstable  behavior  resulted. 

The  degradation  of  performance  due  to  unmodeled  nonlincarities  such  as  actuator  limits  was 
sometimes  found  to  be  worse  than  for  the  nonadaptive  algorithms.  An  adverse  interaction  between 
the  identification  and  the  control  law  developed:  when  the  parameters  generated  by  the  identifi¬ 
cation  began  to  diverge,  control  inputs  actually  caused  the  parameters  to  diverge  more  severely. 
Figure  1.1  illustrates  this  case.  The  left  stabilator  becomes  locked  at  310  seconds.  Steps  of  ±l°/s 
are  made  in  the  pitch  axis  from  305  seconds  to  330  seconds.  This  gives  time  for  the  parameters  to 
adapt  to  the  failure.  At  330  seconds,  the  step  size  is  increased  to  ±5°/s.  In  plots  (a)  and  (b),  the 
adaptive  algorithm  executes  for  the  entire  time  period.  In  plots  (c)  and  (d),  the  adaptive  algorithm 
is  suspended  at  330  seconds,  and  the  control  parameters  remain  constant  beyond  330  seconds.  As 
can  be  seen  in  plot  (b),  when  the  adaptive  algorithm  is  not  suspended,  the  roll  couplings  become 
large  and  even  appear  to  grow  when  the  ±5°/s  step  is  applied.  When  the  adaptive  algorithm  is 
suspended,  the  coupling  in  the  roll  axis  becomes  larger  than  before  the  ±5°/s  step  is  applied  but 
is  smaller  than  the  continuous  adaptive  case  and  remains  within  a  constant  envelope,  as  shown  in 
plot  (d).  Note  that  even  though  the  coupling  is  lower  in  the  suspended  adaptive  case,  increasingly 
poor  tracking  results  during  the  suspended  portion  in  the  pitch  rate  response  in  plot  (c).  Therefore, 
there  is  also  a  trade-off  between  tracking  and  coupling. 

The  improvements  implemented  in  regards  to  instability  are  addressed  in  Chapter  3.  Chapter  4 
proposes  four  methods  of  command  limiting  to  solve  the  actuator  position  and  rate  limiting  problem. 
Chapter  5  addresses  the  issue  of  steady-state  tracking  since  a  steady-state  tracking  error  has  been 
found  to  occur  during  certain  experiments.  Chapter  6  considers  a  more  realistic  flight  simulation 
environment  by  evaluating  the  performance  of  the  aircraft  with  sensor  noise.  The  emphasis  of 
this  thesis  uses  simulations  of  one  flight  condition  and  one  particular  failure.  However,  Chapter  7 
illustrates  the  performance  of  the  adaptive  controller  under  different  flight  conditions  and  failures. 
Finally,  some  miscellaneous  observations  are  made  in  Chapter  8,  and  overall  conclusions  are  made 
in  Chapter  9. 


1.4  Summary 

After  considering  different  types  of  reconfigurable  control,  model  reference  adaptive  control  has  been 
selected  for  the  work  of  this  thesis.  Model  reference  adaptive  control  is  simple,  is  able  to  handle 
cross-couplings,  and  is  of  a  state  feedback  nature.  Previous  work  has  demonstrated  that  model 
reference  adaptive  control  is  useful  for  controlling  an  aircraft  with  a  failure.  However,  robustness 
problems  were  encountered  in  a  broad  range  of  conditions.  The  focus  of  this  thesis  is  to  further 
improve  the  robustness  of  the  adaptive  controller. 
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—  model  response 
- aircraft  response 


Figure  1.1:  Continuous  Versus  Suspended  Adaptive  Algorithm  for  Locked  Stabilator 
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Chapter  2 

Aircraft  Model  and  the  Basic 
Adaptive  Control  Algorithm 


The  purpose  of  this  work  is  to  investigate  reconfigurable  flight  control.  The  aircraft  model  used  in 
the  simulations  and  the  basic  adaptive  control  algorithm  are  discussed  first. 


2.1  Aircraft  Simulation  Model 

The  behavior  of  the  adaptive  controller  has  been  investigated  in  simulations  using  a  model  of  a 
fighter  aircraft  written  in  FORTRAN  by  researchers  at  NASA-Dryden  [7].  The  simulation  uses 
nonlinear  aerodynamic  data  over  the  entire  operational  envelope,  as  well  as  thrust  and  engine 
response  data.  The  aircraft  model  is  that  of  a  high-performance,  supersonic  vehicle  resembling  the 
F-15  fighter  aircraft  shown  in  Figure  2.1,  except  that  the  simulation  assumes  a  single  vertical  tail 
surface  and  a  single  rudder. 

The  aircraft  parameters  are  given  in  Table  2.1.  The  simulation  model  contains  numerous  mod¬ 
ules:  the  equations  of  motion,  the  aerodynamics,  the  atmospheric  model,  the  propulsion  system, 
the  aircraft  actuator  and  surface  command  responses,  the  observation  variables  modeling,  and  a 
basic  proportional-integral  (PI)  control  law  which  may  be  replaced  by  a  user-supplied  control  law. 
The  relationship  between  the  modules  is  shown  in  Figure  2.2. 

The  equations  of  motion  and  the  aircraft  aerodynamics  model  are  six-degree-of-freedom  equa¬ 
tions  representing  the  flight  dynamics  of  a  rigid  aircraft  flying  in  a  stationary  environment  over  a 
flat,  nonrotating  Earth.  The  atmospheric  model  calculates  values  for  speed  of  sound,  acceleration 
due  to  gravity,  air  density,  viscosity,  and  ambient  static  pressure  and  temperature.  The  propulsion 
system  is  comprised  of  two  afterburning  turbofan  engines. 

The  primary  flight  control  surfaces  are  two  horizontal  stabilators,  two  ailerons,  and  a  vertical 
rudder  which  are  shown  in  Figure  2.3.  To  reduce  the  number  of  parameters  that  need  to  be 
identified,  four  commands  were  generated  in  the  original  code:  a  symmetric  command  to  the 
two  stabilators,  a  differential  command  to  the  two  stabilators,  a  differential  command  to  the  two 
ailerons,  and  a  single  command  to  the  one  rudder.  Table  2.2  identifies  the  commands  for  the 
combined  surfaces,  the  limits  imposed  in  the  original  control  law,  and  the  sign  conventions.  The 
combined  commands  are  decomposed  into  individual  surface  commands  given  in  Table  2.3.  Also 
given  are  the  position  and  rate  limits  for  the  commands  to  the  surfaces.  In  this  case,  however,  these 
limits  correspond  to  actual  physical  limits  and  are  implemented  in  the  actuator  simulation  code  as 
opposed  to  the  control  code. 
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Figure  2.1:  F-15  Fighter  Aircraft 


Table  2.1:  Mass  and  Geometry  Characteristics  of  Simulation  Aircraft 


Parameter 

Symbol 

Value 

Wing  Area 

S 

608.0  ft^ 

Wing  Span 

b 

42.8  ft 

Mean  Aerodynamic  Chord 

c 

15.95  ft 

Vehicle  Weight 

Moments  of  Inertia 

w 

45,000.0  lb 

(RoU) 

4 

28, 700.0  slug/ft^ 

(Pitch) 

4 

165, 100.0  slug/ft^ 

(Yaw) 

4 

187, 900.0  slug/ft^ 

Products 

Ixz 

-520.0  slug/ft2 

Ixy 

0.0  slug/ft^ 

lyz 

0.0  slug/ft^ 

Table  2.2:  Command  Input  Definitions  and  Sign  Conventions 


Surface 

Command 

Name 

Symbol 

Supplied 
Control  Law 
Limits  (®) 

Positive 

Sign 

Convention 

Symmetric  Stabilator 

4/c 

+  15/  -25 

trailing  edge  down 

Differential  Stabilator 

Sdc 

±20 

left  trailing  edge  down 

Aileron 

Sac 

±40 

left  trailing  edge  down 

Rudder 

Src 

±30 

trailing  edge  left 

Figure  2.2:  Modular  Structure  of  the  Model 


Table  2.3:  Surface  Deflection  Definition  Equations 


Surface 

Deflection  Definition 

Position  Limit  (®) 

Rate  Limit  (°/s) 

0.55^, 

±20 

24 

^-4  rights 

-0.55,1, 

±20 

24 

•^Sleftc 

0.5(25/f,  +  Sdc) 

±15/-  25 

24 

0.5(25//,  -  6dc) 

±15/-  25 

24 

Src 

^Rc 

±30 

24 

I  y-axis 


Figure  2.3:  Aircraft  Surfaces  and  State  Variables 


Figure  2.4:  Stabilator  Command  Path 

Rate  Position 


Figure  2.5:  Aileron  Command  Path 

The  commands  to  the  individual  actuators  are  generated  from  <5//c,  Sfi^-  The 

command  paths  of  the  stabilators,  ailerons,  and  rudder  are  shown  in  Figures  2.4  ~  2.6,  respectively. 
The  feedback  loops  within  the  figures  are  the  models  for  the  actuators.  The  feedforward  paths  are 
performed  by  the  control  laws.  The  stabilator  path  is  more  complicated  because  of  the  symmetric 
and  differential  capability.  Equivalent  realized  commands  are  also  calculated  based  on  the  actual 
positions  of  each  actuator.  These  equivalent  realized  commands  are  labeled  %,  6d^  and 

The  symmetric  portion  of  the  stabilator  is  primarily  responsible  for  the  pitching  action  in  the 
aircraft.  The  differential  portion  allows  the  stabilator  to  provide  additional  rolling  capability  to 
that  provided  by  the  ailerons.  In  the  original  control  law,  the  command  to  the  differential  stabilator 
was  tied  to  the  aileron  command  with  the  relationship  as  6dc  =  With  this  modification,  the 

control  law  and  identification  are  simplified  by  only  using  three  commands  in  the  control  law:  Snc-i 
Because  of  difficulties  in  dealing  with  actuator  saturation  with  that  choice,  it  was 
decided  to  eliminate  the  coupling  altogether  and  let  5^)^  =  0  for  the  simulations  in  this  thesis  unless 
otherwise  indicated. 

The  observation  variables  are  available  to  the  control  law.  The  supplied  proportional-integral 
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+30  deg/  Rate  Position 

•30  deg  Limit  Limit 

5 

R 

C 


Figure  2.6:  Rudder  Command  Path 

controller  can  adequately  fly  the  aircraft  under  ideal  conditions  but  performs  poorly  in  times  of 
actuator  failure.  The  proportional-integral  controller  was  replaced  with  an  adaptive  control  law 
algorithm. 

2.2  System  Model 

For  the  purpose  of  control  system  design,  a  simplified  model  of  the  aircraft  is  used.  It  is  common 
to  consider  a  reduced-order  model  of  the  form 

X  =  /(x,u), 

y  =  Cx,  (2.1) 

where  x  €  R’'' ,  u  6  R”*,  and  y  €  R^.  The  vector  x  is  chosen  to  include  five  states  of  the  aircraft 
f  1  ^ 

(n  =  5),  with  X  =  a  q  P  p  r  indicating  the  angle  of  attack,  pitch  rate,  angle  of  sideslip, 
roll  rate,  and  yaw  rate,  respectively.  The  interpretation  of  these  variables  is  shown  in  Figure  2.3.  u  is 

[J* 

^Hc  ^Ac  ^Rc  where 

<^.4c)  3.nd  6fic  are  the  commands  to  the  stabilators  (symmetric  portion),  ailerons,  and  rudder, 

r  ^ 

respectively.  It  is  assumed  that  all  the  states  are  measured.  The  output  variable  y  =  q  p  r 
represents  the  variables  to  be  controlled. 

The  trim  state  x*  and  the  trim  control  input  u*  correspond  to  steady,  level  flight  and  are 
obtained  by  solving 

/(x-,u«)  =  0, 

Cx*  =  0.  (2.2) 

A  linearized  model  can  be  derived  with 

X  =  A(x-x*)  +  B(u-u*), 

y  =  Cx.  (2.3) 

With  the  formulation,  x*  and  u*  in  (2.3)  must  be  known  for  control  design.  This  problem  can  be 
alleviated  by  representing  the  system  with  the  following  equations 

X  =  Ax  +  Bu  +  d, 

y  =  Cx,  (2.4) 


11 


where 

d  =  d*  =  -Ax*-Bu*.  (2.5) 

In  this  form,  the  trim  uncertainty  becomes  a  constant  disturbance.  Values  of  the  parameters  A,  B, 
and  d  were  determined  using  a  least-squares  identification  for  four  flight  conditions  and  are  given 
in  Appendix  A.  They  were  also  determined  in  real-time  using  measurements  of  x  and  u  and  for 
control  system  design.  The  trim  values  can  be  determined  from  those  parameters  by  solving 

Ax*-bBu*-fd*  =  0, 

Cx*  =  0.  (2.6) 


In  particular, 


u’  =  -(CB)-^CAx*  -  (CB)"^Cd*, 


(2.7) 


if  the  matrix  CB  is  invertible.  However,  the  control  law  does  not  require  the  explicit  calculation 
of  the  trim  values  to  be  implemented. 


2.3  Control  Law 

The  specific  objective  of  the  control  law  considered  in  this  thesis  is  to  have  the  aircraft  foUow  a 
model  response  specified  by  the  system 

yM=  +  Ba/tw-  (2.8) 

The  reference  input  vm  is  chosen  as  commanded  values  of  angular  rates,  with  tm  =  qc  Pc  Tc  ] 
and  qc,  Pc,  and  rc  represent  rate  commands  in  pitch,  roll,  and  yaw,  respectively.  The  specific 
reference  model  chosen  here  is  of  the  form 

H{s)  =  diag{j^}.  (2.9) 

The  value  of  fc  =  2.5  was  used  in  the  simulations,  following  the  suggestion  of  [8].  Thus,  Aa/  and 
Ba^  in  (2.8)  become 

‘  -2.5  0  0 

Am=  0  -2.5  0  ,  (2.10) 

0  0  -2.5 

■  2.5  0  o’ 

Bm=  0  2.5  0  .  (2.11) 

_  0  0  2.5 

The  control  objective  can  be  achieved  with  the  following  control  law 

u  =  CorAf -b  Gox-b  Vo,  (2.12) 

which  is  shown  in  Figure  2.7.  Indeed,  substituting  (2.12)  in  the  model  given  by  (2.4)  yields 

X  =  Ax  -f  B  •  (CorAf  -b  GqX  -b  Vq)  -b  d 

=  (A -b  BGo)x  "b  BCoTa/ -b  Bvq -b  d.  (2.13) 


12 


Figure  2.7:  Model  Reference  Control  Loop 

Since  y=  C  x, 

y=  (CA  +  CBGo)x  +  (CBCo)rA./  +  CBvq  +  Cd.  (2.14) 

Comparing  the  model  =  AmCxm  +  B\{Tm  with  (2.14)  shows  that  the 

control  objective  can  be  satisfied  with  the  numerical  matrices 

Co  =  (CB)-'Ba/, 

Go  =  (CB)-i(AArC- CA),  (2.15) 

Vo  =  -(CB)-'Cd, 

with  Co  €  Go  G  and  vo  G  For  these  parameters,  the  response  of  the  aircraft 

under  the  control  law  is  identical  to  that  of  the  model.  The  assumption  is  made  that  (CB)"^  is 
finite  [9].  The  transfer  function  of  the  system  (2.4)  is  also  assumed  to  be  minimum  phase  in  order 
to  guarantee  internal  stability  of  the  closed-loop  system. 

2.4  Parameter  Identification  Algorithm 

The  thesis  [6]  proposed  the  use  of  an  input  error  direct  adaptive  control  method.  An  equality, 
proven  in  [6],  was  used  as  the  basis  for  the  identification  technique.  Specifically,  one  has  that 

u  =  CoB^^(y  — A^/y)  +  GqX  +  Vq.  (2.16) 

This  identity  relates  linearly  unknown  parameters  and  known  signals.  Therefore,  standard  least- 
squares  algorithms  can  be  used.  For  example,  a  special  form  of  recursive  least-squares  algorithms 
with  forgetting  factor  was  studied  in  [10]  and  [11],  whose  error  function  is  defined  as 

=  X)  “  0”^[n]w[A:])  ^  A””*-'  +  a  |0[n]  -  e[n  -  1]|^  .  (2.17) 

k=l 
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In  (2.17),  w  and  z  are  vectors  given  by 

(y  -Awy) 

w  =  X  ,  (2.18) 

1 

and 

r  1^ 

z  =  u  =  ^  6hc  ^Ac  ^Rc  \  •  (2.19) 

0  is  a  matrix  of  parameter  estimates  given  by 

r  1 

e=  Gl  .  (2.20) 

L  v?’  J 

The  additional  term  from  the  standard  recursive  least-squares  algorithm  is 
OL  |^[n]  —  6[n  —  l]p.  The  term  penalizes  changes  in  6.  The  least-squares  estimate  corresponding 
to  this  criterion  is  given  by 

=  (52  w[fc]w^[A:]A"“^’  -f  al)~^(^  w[/:]z^[^•]A"“^'  -f-  a6[n  —  1]).  (2.21) 

fc=i  k=l 

In  the  identification,  each  column  of  9  is  separately  calculated  to  determine  the  parameters  used 


in  the  calculations  of  each  of  the  elements  of  the  vector  z.  Defining  P[n]  as 

n 

P[n]  =  w[A:]w^[A:]A"-^-  +  al)"' ,  (2.22) 

k-l 

a  recursive  form  for  the  inverse  of  the  covariance  matrix  P[7i]  is 

P“i[n]  =  AP“^[n  -  1]  -f  w[n]w^[7i]  -|-  a(l  -  A)I,  (2.23) 

with  the  initial  condition  that  P~^[0]  =  al.  The  recursive  formula  for  0[n]  is 

0[n]  =  9[n  -  1]  -f  P[n]w[n](z^[n]  -  w^[n]0[n  -  1])  +  aAP[n](0[Ti  -  1]  -  9[n  -  2]).  (2.24) 

A  recursive  form  of  P[n]  is  also  desirable,  so  that  P“^[n]  does  not  need  to  be  inverted.  An 
approximate  recursive  form  of  P[n]  was  proposed  in  [11]  as 

PN  =  ;^(P[^  -  1]  -  P(«  -  l]C[n](AI  -h  C^[n]P[n  -  l]C[n])-'C^[n]P[n  -  1]),  (2.25) 


with  C[n]  =  w[n]  \/ma(l  —  A)e,-  |,  where  m  is  the  dimension  of  w[ti],  and  e,-  is  a  vector  of 

zeros  except  the  position  which  is  one.  The  index  i  is  incremented  at  each  iteration  and  returns 
to  one  when  the  end  of  the  vector  is  reached.  A  well-known  problem  in  adaptive  control  is  that 
precision  errors  in  calculating  P[n]  can  result  in  P[n]  not  necessarily  being  symmetric.  Therefore, 
the  symmetry  of  P[n]  was  used  and  enforced  at  each  iteration.  The  stabilized  recursive  least-squares 
algorithm  given  by  (2.24)  and  (2.25)  was  found  to  alleviate  problems  occurring  during  periods  of 
low  excitation,  and  the  appro.ximate  form  (2.25)  was  found  to  perform  satisfactorily. 
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2.5  Performance  Evaluation 


Two  types  of  failures  were  simulated  in  the  code.  The  first  failure  is  a  locked  surface  such  that 
the  surface  no  longer  responds  to  the  commands  and  stays  fixed.  The  second  failure  is  a  missing 
stabilator  surface.  A  missing  stabilator  is  a  more  severe  failure  than  a  locked  stabilator  because 
not  only  is  control  power  reduced,  but  also  the  stabilization  effect  of  the  tail  surface  is  lost.  In  the 
code,  a  missing  stabilator  was  simulated  by  having  the  surface  follow  the  direction  of  the  local  angle 
of  attack.  The  local  angle  of  attack  is  not  exactly  known,  but  was  approximately  calculated  as  in 
Appendix  B.  Both  the  locked  and  missing  failures  were  separately  applied  to  the  left  stabilator 
in  the  simulations.  A  locked  aileron  failure  was  also  simulated  independently,  and  together  with  a 
locked  stabilator. 

The  main  flight  condition  for  the  simulations  is  an  altitude  of  9,800  feet  and  a  speed  of  Mach  0.5. 
In  addition  to  the  main  flight  condition,  three  additional  flight  conditions  are  considered  in  Chap¬ 
ter  7,  as  suggested  by  [7]  to  cover  the  operational  envelope: 

•  9,800  feet,  Mach  0.9, 

•  39,800  feet,  Mach  0.6,  and 

•  39,800  feet,  Mach  1.4. 

The  proportional-integral  controller  provided  with  the  original  code  is  able  to  maintain  flight  for 
the  unfailed  aircraft.  In  the  case  of  the  failed  aircraft,  the  proportional-integral  controller  only 
works  over  limited  periods  of  time.  In  both  cases,  however,  responses  of  the  aircraft  to  brief 
maneuvers  can  be  used  to  determine  the  parameters  of  a  linearized  model  using  batch  least-squares 
algorithms.  Two  different  sets  of  the  parameter  matrices  Cq,  Gq,  and  vq,  corresponding  to  the 
unfailed  and  locked  stabilator  conditions,  were  generated  using  a  10  second  data  sample.  To  provide 
excitation,  steps  were  applied  to  the  proportional-integral  controller  first  in  altitude,  next  in  roll 
angle,  and  finally  in  the  angle  of  sideslip.  The  details  of  the  identification  procedure  are  found 
in  Appendix  A.  The  results  of  these  identification  experiments  are  useful  to  evaluate  the  baseline 
controller  performance  independent  of  the  adaptation  response,  as  well  as  to  determine  nominal 
controller  parameter  values.  However,  the  recursive  algorithms  do  not  use  the  off-line  identification 
results  for  the  failed  aircraft.  The  state-space  matrices  are  given  for  the  four  flight  conditions 
in  Appendix  A  for  both  the  unfailed  and  locked  stabilator  aircraft  with  the  differential  stabilator 
movement  enabled  (Sdc  =  Appendix  A  also  contains  the  state-space  matrices  for  the  unfailed 

and  locked  stabilator  aircraft  with  the  differential  stabilator  movement  disabled  (Sp^  =  0). 

The  performance  of  the  control  law  is  evaluated  throughout  the  thesis  by  applying  steps  of 
commands  in  pitch  rate  or  roll  rate.  The  control  law  is  activated  at  305  seconds  for  both  the 
batch  and  recursive  identifications.  The  selection  of  305  seconds  is  arbitrary,  but  a  period  of  time 
is  necessary  to  ensure  that  all  transients  of  the  simulation  code  have  decayed.  Failures  occur  at 
305  seconds  for  the  batch  experiments  and  310  seconds  for  the  adaptive  experiments.  Runs  are 
typically  evaluated  for  50  seconds.  The  aircraft  code  calculates  the  outputs  of  the  aircraft  at  a 
rate  of  50  Hz.  The  same  sampling  rate  is  used  in  the  identification  and  the  control  law.  For 
the  recursive  algorithms,  the  parameter  matrices  are  initialized  with  the  values  from  the  unfailed 
identification  case.  Design  parameters  of  the  RLS  controller  are  typically  set  to  A  =  0.99  and 
a  =  10.0.  Appendix  C  provides  the  details  of  all  the  simulations  in  the  thesis. 
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2.6  Summary 


This  chapter  presented  the  aircraft  simulation  model  and  the  adaptive  control  algorithm  that  were 
used  in  the  thesis,  based  on  previous  work.  A  reduced-order  model  simplified  the  highly-detailed, 
nonlinear  model  of  a  fighter  aircraft.  The  adaptive  algorithm  used  a  model  reference  adaptive 
control  law  which  commanded  the  aircraft  to  follow  a  first-order  response  in  angular  rates.  The 
model  was  further  simplified  by  designing  the  control  law  to  calculate  only  three  common  control 
inputs.  In  this  way,  the  number  of  parameters  that  needed  to  be  identified  was  reduced.  Parameters 
were  identified  using  an  indirect  error  equation  technique  and  a  stabilized  recursive  least-squares 
procedure.  The  performance  of  the  adaptive  control  algorithm  was  demonstrated  through  the 
aircraft  simulations. 
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Chapter  3 

Modifications  to  Improve  Robustness 


The  formulation  of  the  control  law  is  given  by  (2.12).  The  matrices  Co,  Go,  and  vq  are  replaced  by 
the  estimates  provided  by  the  stabilized  recursive  least-squares  algorithm  with  forgetting  factor.  In 
the  course  of  the  research,  certain  problems  developed  with  this  implementation.  New  formulations 
to  identify  parameters  were  created  and  found  to  eliminate  many  of  these  problems.  First,  an 
inverse  controller  formulation  was  developed  to  address  a  singularity  problem  occurring  in  the 
Co  matrix.  In  some  cases,  the  Co  matrix  approached  singularity,  and  instability  occurred  shortly 
thereafter.  The  problem  was  largely  eliminated  by  identifying  the  inverse  of  the  Co  matrix,  instead 
of  Co.  Eventually,  this  formulation  was  transformed  into  a  form  of  indirect  identification  which  led 
to  the  best  results  overall. 

Other  modifications  include  the  filtering  of  the  actuator  commands,  to  account  for  the  delay  in 
the  actuator  responses,  and  the  normalization  in  the  algorithm,  to  improve  responses  with  demand¬ 
ing  flight  maneuvers.  Another  modification,  the  use  of  actuator  measurements  in  the  identification, 
was  experimented  with  but  not  retained. 


3.1  Indirect  Identification  Formulation 


The  adaptive  algorithm  is  not  guaranteed  to  be  stable  unless  certain  modifications  are  incorporated 
to  ensure  that  the  Cq  matrix  is  nonsingular  [12].  Unfortunately,  these  modifications  are  very  difficult 
to  implement.  A  simple  modification  that  was  found  to  alleviate  the  problems  is  that  of  an  “inverse” 
identification  formulation.  The  inverse  Cj  of  the  Cq  matrix  is  identified.  The  estimate  matrix  Ci 
is  then  inverted  during  each  iteration  of  the  algorithm  to  give  Cq  which  is  used  in  the  control  law. 
The  original  error  equation  is  of  the  form 


Wo  =  CJWi  +  GqV/2  +  Vq, 


where 


Wo 


%c  ^Ac  ^Rc 


'vi  =  B^/(y-AMy), 

W2  =  \  a  q  /3  p  r 


(3.1) 


(3.2) 


A  new  equation  is  generated  with  the  output  based  on  the  model  and  the  surface  commands  being 
part  of  the  input.  Essentially,  wo  and  wi  are  “swapping”  positions  in  (3.1).  The  new  equation  is 
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defined  as 


wi  =  Cj-wo  +  G1W2  +  Vj.  (3.3) 

In  order  to  obtain  the  controller  parameters  from  the  adaptive  parameters,  the  relationship  between 

the  nominal  parameters  is  used: 

Co  =  Cj  \ 

Go  =  -C^^Gi,  (3-4) 

Vo  = 

The  resulting  values  for  Co,  Go,  and  vo  are  applied  to  the  control  law  (2.12).  Using  the  original 
formulation,  the  vector  w[n]  and  z[n]  are  defined  by  (2.18)  and  (2.19),  respectively.  Using  the 
inverse  formulation,  one  has  instead 

'  %  ' 

w[n]  =  ,  (3.5) 

X 
1 

and 

z[7i]  =  Bj}{y  -Amy).  (3.6) 

The  inverse  formulation  is  more  complex  than  the  original  formulation  primarily  because 
must  be  calculated  at  each  iteration.  Further,  Ci  could  be  singular  in  this  algorithm,  so  that 
Co  would  be  unbounded.  However,  this  situation  was  not  found  to  occur  in  simulations,  so  that 
the  additional  computations  were  worth  the  cost.  Simulations  showed  that  the  inverse  control  law 
formulation  performs  better  than  the  original  control  law  formulation  over  various  conditions. 

The  inverse  controller  formulation  is  closely  related  to  a  special  form  of  indirect  identification 
applicable  to  the  model  reference  control  problem.  Typically,  for  an  indirect  identification,  the 
state-space  model 

x=Ax  +  Bu  +  d  (3.7) 

is  used  to  estimate  the  parameters  A,  B,  and  d.  However,  using  (3.7)  requires  the  identification 
of  (n  +  m  -f  1)  X  n  parameters,  which  is  substantially  larger  than  that  required  by  the  previous 
methods.  Instead  of  using  (3.7)  in  the  identification, 

y=  C  x=  CAx  +  CBu  +  Cd  (3.8) 

can  be  used  to  determine  CA,  CB,  and  Cd,  which  are  the  matrices  needed  to  calculate  the 
controller  parameters,  and  then  only  (n  +  m  +  1)  X  m  parameters  are  required.  For  n  =  5  and 
7n  =  3,  45  parameters  are  required  using  (3.7),  and  27  parameters  are  required  using  (3.8).  In  the 
latter  case,  the  number  of  parameters  is  the  same  as  the  number  of  controller  parameters.  From 
the  matrices  CA,  CB,  and  Cd,  the  matrices  Co,  Go,  and  vo  can  be  calculated  as  in  (2.15),  that  is 

Co  =  (CB)-IBa/, 

Go  =  (CB)-1(Aa/C-CA),  (3.9) 

Vo  =  -(CB)"^Cd. 
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For  this  identification  formulation,  the  vectors  required  by  the  identification  algorithm  are 


w[n]  = 


^Hc 

^Ac 

SRc 

X 

1 


(3.10) 


and 


z[n]  =y= 


T 


q  p  r 


(3.11) 


Note  that  the  identification  procedure  no  longer  depends  on  the  model  parameters  Am  and  Bm- 
This  fact  is  useful  for  the  development  of  the  command  limiting  routines  presented  in  Chapter  4. 
Otherwise,  the  results  obtained  with  this  method  were  found  to  be  identical  to  those  obtained  with 
the  inverse  controller  formulation. 


3.2  Filtering  of  Actuator  Commands 

Delays  within  the  system  can  cause  instabilities.  One  such  location  for  delays  is  in  the  actuator 
command  path.  As  shown  in  Figures  2.4  -  2.6,  the  equivalent  realized  commands  (the  outputs  of  the 
systems)  and  the  surface  commands  (the  inputs  to  the  systems)  are  not  identical,  even  if  saturation 
is  not  encountered.  The  unity  feedback  system  with  an  open  loop  containing  an  integrator  and 
limiting  routines  modifies  the  surface  command.  If  the  position  and  rate  limits  are  not  encountered, 
the  command  path  is  modeled  in  the  aircraft  code  by  the  first-order  system 

(3.12) 

To  increase  the  accuracy  of  the  identification,  the  control  variables  (5//^,  S^c,  and  in  (3.5)  or 
(3.10)  were  filtered  in  the  control  code  using  (3.12).  A  zero-order  hold  approximation  was  used, 
as  opposed  to  the  Euler  approximation  implemented  in  the  actuator  code,  so  that  the  adaptive 
algorithm  did  not  have  available  an  exact  replica  of  the  actuator  dynamics.  Figure  3.1  illustrates 
that  a  low  condition  number  of  the  CB  matrix  and  stability  result  with  the  inclusion  of  the  surface 
command  filtering.  Otherwise,  instability  rapidly  occurs.  The  filtered  values  of  the  commands  are 
only  applied  to  the  identification  procedure.  The  unfiltered  values  of  Shc,  and  are  passed 
as  inputs  to  the  systems  shown  in  Figures  2.4  -  2.6. 


3.3  Use  of  Actuator  Measurements 

Another  modification  that  contributes  to  the  solution  of  the  instability  problem  is  utilizing  the 
actuator  measurements  in  place  of  the  commands  to  the  actuators.  In  other  words,  instead  of 
using  the  command  control  variables  6hc^  Sac,  and  in  (3.5)  or  (3.10),  the  equivalent  realized 
commands  Sff,  6a,  and  6^  are  utilized.  Figure  3.2  demonstrates  the  effect  of  using  actuator  mea¬ 
surements  versus  using  the  commands  to  the  actuators  for  a  missing  stabilator  failure.  Results  are 
similar  between  the  filtering  of  actuator  commands  and  the  use  actuator  measurements  techniques. 
Although  it  would  seem  preferable  to  use  the  actuator  measurements  if  they  are  available,  the 
procedure  has  a  drawback.  Any  change  in  the  actuator  path,  such  as  a  locking  of  a  stabilator,  will 
not  be  identified  by  the  procedure  and  will  have  to  be  separately  detected.  With  actuator  filtering, 
the  response  of  the  actuators,  surfaces,  and  aircraft  are  lumped  together  and  the  overall  system 
dynamics  are  identified. 
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Pitch  Rato  Response  to  PKch  Rale  Command:  UnfiRered  &  FiRered  Commands  Condition  of  C8  Matrix:  UnfiHered  &  FiRered  Commands 


—  model  response 

- unfiltered  commands 

—  •  —  filtered  commands 

Figure  3.1:  Unfiltered  versus  Filtered  Commands  for  Pitch  Rate  Response  to  Pitch  Rate  Command 


3.4  Normalization  of  the  Algorithm 

Instability  worsens  as  the  demand  of  the  required  task  increases.  For  example,  steps  of  ±l°/s  in 
many  cases  did  not  generate  unstable  responses  in  the  simulations.  However,  steps  of  ±5®/s  or 
±10^/s  often  resulted  in  instability.  A  contributing  factor  to  this  instability  appears  to  be  that, 
as  the  numbers  in  the  control  algorithm  become  large,  the  effect  of  the  unmodeled  dynamics  is 
increased.  A  solution  to  this  problem  is  to  include  normalization  in  the  algorithm.  The  vectors  w[n] 
and  z[n]  defined  by  (3.5)  and  (3.6),  respectively,  or  (3.10)  and  (3.11),  respectively,  in  the  recursive 

least-squares  algorithm  are  normalized  by  dividing  each  element  of  the  vectors  by  ^w^[n]w[n]. 
/?  is  a  user-defined  constant.  Simulations  have  shown  that  ^  =  0.1  yielded  good  results  for  the 
system. 

Figure  3.3  illustrates  the  effect  of  including  normalization  in  the  algorithm  for  a  missing  sta- 
bilator  failure.  The  glitches  in  the  pitch  rate  command  response  in  plot  (a)  are  reduced  using  the 
normalization.  The  spiked  responses  in  the  roll  rate  coupling  shown  in  plot  (b)  are  also  reduced. 
The  normalization,  however,  is  not  found  by  itself  to  resolve  all  of  the  problems  associated  with 
large  maneuvers  and  due  to  actuator  position  and  rate  limits.  More  complete  solutions  are  discussed 
in  the  next  chapter. 

3.5  Summary 

Instability  problems  appeared  in  the  course  of  the  research.  Simulations  showed  that  the  method 
by  which  the  parameters  were  identified  affected  the  occurrence  of  instability.  Two  alternate  iden¬ 
tification  formulations  were  proposed  to  reduce  the  likelihood  of  instability.  The  inverse  controller 
formulation  is  a  direct  approach  that  uses  the  same  control  law  as  the  original  identification  but 
requires  the  calculation  of  a  matrix  inverse  to  obtain  the  parameter  matrices.  The  indirect  iden- 
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Right  Stabilator  Position  (deg)  q  (dog/t) 


Pitch  Rsta  Rasponsa  to  Pitch  Rate  Command:  Using  Actuator  Maasuramants  Roll  Rate  Response  to  Pitch  Rate  Command;  Using  Actuator  Maasuraments 


Right  Stabilator  Position  lor  Pitch  Rale  Command;  Using  Actuator  Measurements  Right  Aileron  Position  for  Pitch  Rata  Command:  Using  Actuator  Maasuramants 


Tima  (sac)  Time  (sac) 

(c)  (d) 

—  model  response 

- using  filtered  commands 

—  '  —  using  actuator  measurements 

Figure  3.2:  Results  Using  Actuator  Measurements  for  Missing  Stabilator  Failure 
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tification  formulation  is  closely  related,  with  the  added  benefit  of  removing  the  model  parameters 
from  the  identification.  Instability  was  also  eliminated  through  the  use  of  filtering  of  the  actuator 
commands,  to  more  accurately  account  for  delays  within  the  actuators,  and  normalization  of  the 
algorithm,  to  reduce  the  dynamic  range  of  the  signals.  The  use  of  actuator  measurements  improves 
the  robustness  of  the  system  but  has  a  drawback  because  it  requires  additional  fault  detection  logic 
to  handle  failures  related  to  the  actuators. 
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Chapter  4 

Actuator  Position  and  Rate 
Saturation 


The  performance  of  the  adaptive  control  law  discussed  so  far  deteriorates  rapidly  when  the  size  of 
the  steps  of  the  reference  input  increases  from  ±F/s  to  ±10®/s.  Figures  4.1  and  4.2  illustrate  the 
response  of  the  aircraft  to  pitch  rate  commands  of  step  sizes  ±l®/s  and  ±10°/s,  respectively.  In 
the  large  step  size  responses,  the  response  of  the  aircraft  due  to  the  pitch  rate  command  does  not 
match  the  model  and  has  glitches  that  do  not  appear  in  the  small  step  size  plots.  Of  troubling 
concern  is  the  large  coupling  up  to  ±10°/s  in  the  roll  axis,  shown  in  Figure  4.2  (b),  which  does 
not  decay  with  time.  Part  of  the  problem  is  that  the  actuators  controlling  different  axes  do  not 
saturate  in  a  coherent  way,  so  that  decoupling  is  lost.  Therefore,  the  problem  is  not  necessarily  one 
of  insufficient  control  power  but  of  proper  usage  of  the  control  power  available.  Rate  saturation 
may  be  observed  in  the  right  stabilator  and  aileron  responses.  Plots  (c)  and  (d)  of  Figure  4.2  show 
the  right  stabilator  and  aileron  positions,  respectively.  Linear  portions  in  the  position  plots  are  due 
to  extensive  rate  limiting  to  ±24°/s.  Plots  (a)  and  (b)  of  Figure  4.3  show  the  corresponding  rates. 

There  are  two  options  to  handle  the  saturation  problem.  The  first  is  to  design  a  control  law 
that  takes  the  constraints  into  account.  For  example,  a  constrained-optimization  procedure  may  be 
applied  to  calculate  the  values  of  the  control  inputs  needed  to  optimize  some  cost  criterion  subject 
to  the  constraints.  The  second  is  to  develop  a  control  law  without  considering  the  constraints  in 
the  design  of  the  control  law  and  then  devising  a  procedure  to  modify  the  control  signals  when 
the  constraints  are  not  satisfied.  For  example,  admissible  control  inputs  may  be  calculated  to 
approximate  the  desired  control  inputs  in  some  optimal  sense  within  the  constraints.  The  first 
option  is  more  difficult  to  develop  and  requires  complex  calculations.  The  second  option,  generally 
known  as  command  limiting,  has  the  advantages  that  the  nominal  control  system  design  is  separate 
from  the  treatment  of  control  saturation  and  that  the  resulting  algorithm  is  of  low  to  moderate 
complexity.  In  general,  command  limiting  is  related  to  control  allocation,  which  is  the  problem 
of  distributing  control  requirements  among  multiple,  constrained  actuators.  Control  allocation 
has  been  studied  in  [13]  -  [15],  whose  concept  is  to  find  an  allowable  control  input  such  that  the 
generated  moment  in  the  direction  of  the  desired  moment  is  closest  to  the  desired  value.  In  [8], 
the  idea  is  to  find  an  approximation  of  the  acceleration  that  is  closest  to  the  desired  value  in  a 
least-squares  sense.  In  [16],  a  method  is  proposed  where  the  reference  input  is  approximated  so 
that  the  modified  reference  input  applied  to  the  control  law  does  not  yield  saturation.  Four  control 
allocation  methods  are  proposed  in  this  thesis  based  on  these  concepts  to  address  the  saturation 
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Pitch  Ra(«  Responae  to  Pitch  Rato  Command:  Modot  &  No  Command  Limiting 


(a) 


RigN  Stabitalor  Position  for  Pitch  Rato  Command:  No  Command  Limiting 


(c) 


(b) 


(d) 


model  response 
aircraft  response 


Figure  4.1:  Responses  to  a  Small  Pitch  Rate  Command  without  Command  Limiting 
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Right  SUbilalor  Position  (deg)  q  (deg/s) 


Pitch  Rale  Response  to  Pitch  Rate  Command:  Model  &  No  Command  Limiting 


Roll  Rate  Response  to  Pitch  Rate  Command:  Model  &  No  Command  Limiting 


—  model  response 
- aircraft  response 


Figure  4.2:  Responses  to  a  Large  Pitch  Rate  Command  without  Command  Limiting 
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model  response 
aircraft  response 


Figure  4.3:  Actuator  Rates  for  a  Large  Pitch  Rate  Command 


problem.  These  methods  are 

•  scaling  of  control  inputs, 

•  relaxation  of  control  requirements, 

•  scaling  of  reference  inputs,  and 

•  a  least-squares  approximation  to  the  commanded  accelerations. 

The  methods  are  compared  by  considering  computations,  ease  of  implementation,  and  performance 
as  measured  by  simulations.  The  presentation  follows  the  organization  of  the  paper  by  Bodson  and 
Pohlchuck  [17].  Although  the  methods  were  proposed  earlier,  their  evaluation  in  a  reconfigurable 
control  application  is  novel. 


4.1  Definition  of  Constraints 
r  1^ 

The  input  u  =  8hc  ^Ac  ^Rc  must  be  in  the  admissible  position  control  set  Up  where  Up  = 
{u  I  for  z  =  1 . .  .3,Pt^Tnm  ^  <  Pi, max}-  The  rate  of  variation  of  u  is  also  constrained  so  that 

I  I  <  di^rnax  ioT  i  =  1...3.  Because  the  application  of  digital  flight  control  is  considered,  the 
constraint  on  the  rate  of  variation  can  be  translated  into  a  position  constraint.  Given  a  control  input 
u(n  -  1)  at  the  previous  time  sample  and  a  sampling  period  T,  the  admissible  control  rate  set  is 
defined  to  be  Ur,  where  =  {u  |  for  z  =  1 . .  .3,  u,(n-  l)  +  Tdi^max}- 

Both  the  position  and  rate  limits  can  be  translated  into  a  single  admissible  control  set,  U  =  UpOUr, 
with  U  of  the  form  =  {u  [  for  Z  =  1  .  .  ^S^Ui^rnin  <  '^i  ^  ^t,Tnax}- 


27 


4.2  Method  #1:  Scaling  of  Control  Inputs 

The  first  method  is  the  simplest.  The  desired  control  input  is  denoted  and  the  control  input 
produced  by  the  command  limiting  technique  is  u.  The  idea  behind  this  method  is  to  scale  back 
the  commands  proportionally  so  that  the  constraints  are  satisfied.  First,  position  constraints  are 
checked.  Of  the  three  control  elements,  the  one  which  exceeds  a  limit  by  the  largest  amount  (if 
any)  determines  the  scaling  factor  which  will  proportionally  reduce  the  values  of  the  other  two 
elements.  All  the  control  elements  will  belong  to  the  admissible  position  set  after  the  scaling  factor 
is  applied  to  all  three  commands.  Second,  the  rate  constraints  are  checked,  already  taking  into 
account  the  position  limiting  adjustment.  All  three  control  elements  are  checked  to  see  if  a  rate 
limit  is  exceeded.  A  scaling  factor  is  determined  which  reduces  the  largest  rate  excess  to  within 
acceptable  limits.  The  same  scaling  factor  is  applied  to  the  other  two  control  elements  so  that  all 
three  elements  belong  to  the  admissible  rate  set.  Note  that  by  convexity  of  the  admissible  position 
set,  it  follows  that  the  final  control  value  still  belongs  to  the  admissible  position  set.  The  algorithm 
is  described  below. 

Algorithm  #1: 

1.  Let  Pi  be  the  largest  number  such  that  0  <  pi  <  1  and  ui(7i)  =  piUd{n)  G  Up. 

2.  Let  p2  be  the  largest  number  such  that  0  <  p2  <  1  and  U2(n)  =  u(n  —  1)  +  P2(piUrf(n)  — 
u(n  “  1))  G  Ur.  Let  u(n)  =  U2(n). 

The  position  limits  are  evaluated  before  the  rate  limits.  The  idea  is  that  performing  the  limiting 
in  this  order  is  more  likely  to  preserve  the  directionality  of  the  commands. 

4.3  Method  ^2:  Relaxation  of  Control  Requirements 

The  second  method  allows  the  control  requirements  to  be  relaxed  when  the  constraints  are  violated. 
To  implement  this  method,  an  alternate  expression  is  used  for  the  control  law.  The  original  control 
law  is 

Ud  =  CqVm  +  Gqx  +  Vo  (4.1) 

Substituting  the  values  for  Cq,  Gq,  and  vq  as  given  by  (2.15)  in  (4.1)  yields  the  control  law 

Ud  =  (CB)-^BA/rM  +  (CB)-*(AMC-CA)x-(CB)-'Cd 

=  (CB)-^BA/rM  +  (CB)-^AA.fCx-  (CB)-’CAx-  (CB)-^Cd 

=  (CB)-^BwrM  +  (CB)-iAA/y-(CB)-iCAx-(CB)-^Cd  (4.2) 

By  setting  —Am  =  —  kl,  (4.2)  becomes 

Ud  =  (CB)-U’riv/-(CB)-U’y-(CB)-'CAx-(CB)-‘Cd 

=  (CB)-U-(rjv/-y)  +  (CB)-^(-CAx-Cd)  (4.3) 

In  (4.3),  Ud  depends  linearly  on  fc.  The  control  requirements  are  primarily  embedded  in  the  con> 
stant  k  and  can  be  relaxed  by  reducing  its  value.  Let  the  control  vector  be  decomposed  as  the 
sum  of  two  vectors  and  u^,  given  by 


=  Ua  +  Ufc, 

Ua 

=  {CB)~^k{TM-y), 

(4.4) 

U6 

=  (CB)-i(-CAx- Cd). 
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(4.4)  gives  a  partition  of  (4.3).  The  idea  of  the  method  is  to  find  a  p  that  satisfies  0  <  p  <  1  and 
puo  +  Ube  U.  In  essence,  p  is  reducing  the  value  of  k  in  the  equation  for  Ua  in  (4.4).  This  method 
is  similar  to  Algorithm  but  the  coding  is  complicated  by  the  term. 

Algorithm 


1.  Find  the  range  Ii,\  such  that  ^  “b  ^6,1  ^  for  ull  p  ^ 

[ri.mtnri.TOor]-  The  range  can  be  obtained  by  considering  9  cases: 


(a) 


(b) 

(c) 

(d) 

(e) 

(f) 

(g) 

(h) 

0) 


if  '^b,l  ^  Hi, min  U.nd  Uap  *b  Ub,l  ^  '^l,maxj  lot 


_  ^6,1  _  ~  '^^6,1 
PlfTTitn  —  )  Pi, mar  —  ? 

'^a,l  ^a,l 


(4.5) 


and  fi^min  —  Pi, mini  1*1, mas  —  Pi, maxi 

if  lib,!  ^  ^l,min  3'Ild  Ux^fnin  ^  ^a,l  "f"  “^b,!  ^  ^l,maxi  lot  ^l,max  —  1  and  Vx^min  ~  Pl,mini 

where  pi.min  is  defined  in  (4.5); 

if  Ubp  <  111, min  and  Ua.i  +  M6,i  <  ^‘■i,mini  Ri  =  0,  the  empty  set; 

if  Hi, min  <  H6,1  <  Hi, mar  and  Ua,i  +  U(,,i  >  Ux,maxi  let  ri,min  =  0,  ri,max  =  Pl,maxi  where 

Pl.mos  is  defined  in  (4.5); 

if  Hi,min  ^  Hj,i  <  ll‘l,max  and  Ux^rnin  ^  Ho,i  +  U;,,!  ^  'f^l,maxi  let  Vx^min  —  Oj  I'l.mas  =  1; 

if  Hi,min  <  H(,,i  <  Ux,max  and  Ua,l  +  H6,l  <  Hi,min)  let  ri,mm  =  0,  ri,mai  =  Pi, mini  where 

Pi, min  is  defined  in  (4.5); 

if  lit,!  >  ui,max  and  Ua,i  +  Ub,x  >  Hi,max,  Ri  =  0,  the  empty  set; 

if  Ub,x  >  Ux,max  and  Ul,min  <  «„,!  +  Ubp  <  Hi,mai,  let  ri,max  =  1,  ri,min  =  Pl,max,  where 

Pi,mox  is  defined  in  (4.5); 

if  U6,1  >  Hi,max  and  Ua,l  +  H6,i  <  Ul,min,  let  ri,min  =  Pi, max,  l’i,mai  =  Pl,min,  where 
Pi,min,  Pi,max  are  defined  in  (4.5). 


2.  Repeat  step  1  with  appropriate  index  changes  to  determine  the  ranges  of  R2  and  R^. 

3.  Find  R  =  RxO R2(^ R^-  If  R  is  not  empty,  let  p  be  the  largest  number  in  R  and  u  =  pua  +  U(,. 
If  R  is  empty,  let  u  be  such  that,  for  i  =  1 . .  .3,  u,-  =  w,-,min  if  hj,,-  <  Ui,min,  h,-  =  Ui^max  if 
H(i,i  >  Ui^maxi  aiid  U{  =  Ud,i  Otherwise. 


The  tests  are  based  on  the  values  of  Ub,x  and  Ua,i  +  H6,i  which  are  the  values  of  pua,x  +  H6,i  for  p  =  0 
and  p  =  1,  respectively.  The  variables  pi,min  and  px,max  are  the  values  of  p  such  that  pua,i  +  H6,i 
is  equal  to  Ux,min  and  Hi,max,  respectively.  These  values  may  be  outside  the  allowable  interval,  but 
are  not  calculated  in  that  case.  Division  by  zero  in  (4.5)  is  avoided. 

A  scaling  factor  cannot  be  determined  in  the  case  of  R  being  empty  in  step  3.  This  condition  is 
handled  by  applying  a  simple  saturation  function.  The  case  occurred  occasionally  in  simulations, 
but  not  often. 


4.4  Method  #3:  Scaling  of  Reference  Inputs 

This  method  is  similar  to  the  second  method,  in  that  an  input  signal  is  scaled.  Instead  of  scaling 
the  control  vector  u^,  the  reference  input  rA/(n)  is  scaled.  This  method  is  computationally  close  to 
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the  second.  Indeed,  Uo  and  u;,  were  defined  in  (4.4).  and  uj,  can  be  redefined  by  the  following 
equations: 


U<i  =  Ua  +  U6, 

Ua  =  (CB)"UrA/,  (4.6) 

U6  =  (CB)-i(-CAx- Cd-fty). 

The  same  procedure  used  in  Method  ^2  is  applied  to  this  method.  Again,  the  control  vector  u  is 
determined  such  that  0  <  /?  <  1  and  pUa  +  £  U.  A  difficulty  with  this  technique  is  that  it  is 

far  more  common  to  be  in  a  situation  where  no  value  of  p  satisfies  the  requirements.  This  problem 
particularly  occurs  when  the  reference  input  moves  from  a  large  value  to  a  zero  value.  In  this  case, 
scaling  does  not  make  any  difference.  To  resolve  the  problem,  the  algorithm  can  be  modified  so 
that  a  linear  combination  of  the  previous  reference  input  and  the  current  reference  input  is  used. 
Instead  of  replacing  rA./(7i)  by  /)riv/(n)  in  given  by  (4.6),  rA/(n)  is  replaced  by  the  following 
equation: 

rM{n)  =  pVMin)  +  (1  -  p)rM{n  -  !)■  (4.7) 

For  /)  =  1,  the  modified  reference  input  is  equal  to  the  current  reference  input.  For  p  =  0,  the 
modified  reference  input  is  equal  to  the  previous  reference  input.  Algorithm  ^2  can  again  be 
applied,  but  with  changes  to  the  definitions  of  and  uj  in  (4.6)  given  in  the  following  equation: 

uj  =  Ua  +  U6, 

Ua  =  (CB)-UTA/(n)  -  (CB)-U'rA/(n  -  1),  (4.8) 

U6  =  (CB)“^(-CAx- Cd- A:y)  + (CB)-UTAi(n  -  1). 

The  term  (CB)“^A:rA/(n—  1)  is  subtracted  from  Ua  but  is  also  added  to  u;,.  The  resulting  equation 
for  uj  is  the  same  as  in  (4.6).  The  algorithm  is  described  below. 

Algorithm  ^3: 

1.  Find  the  largest  p  such  that  0  <  p  <  1  and  pUa  +  uj  G  17,  using  the  method  of  Algorithm  #2, 
but  with  Ua  and  u;,  defined  in  (4.6).  If  a  p  is  found,  let  u  =  pUa  +  uj  and  stop.  Else,  proceed. 

2.  Find  the  largest  p  such  that  0  <  p  <  1  and  pUa  +  uj  G  17,  using  the  method  of  Algorithm  ^2, 

but  with  Ua  and  uj  defined  in  (4.8).  If  a  p  is  found,  let  u  =  pUa  +  uj,  and  replace  rA^f(n)  using 
(4.7).  If  no  p  can  be  found,  let  u  be  such  that,  for  i  =  1 . .  .3,  u,-  =  if  uj.i  <  Uj.mtn) 

n-i  =  Ui.max  if  >  y-i,max,  and  u;  =  Ud,i  Otherwise. 

Cases  can  still  arise  where  a  value  of  p  cannot  be  determined.  If  the  states  and  parameters 
vary  by  large  amounts  during  a  sampling  period,  scaling  of  the  reference  input  may  not  bring  the 
control  input  within  the  constraint.  The  simple  saturation  function  described  in  Method  ^2  is 
again  applied. 

4,5  Method  ^4:  LS  Approximation  of  Commanded  Accelera¬ 
tions 

This  method  determines  the  command  limiting  based  on  the  commanded  accelerations.  The  accel¬ 
eration  produced  is  determined  from  the  state-space  model 

y=CAx-hCBu-l-CBd,  (4.9) 
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and  the  component  due  to  the  desired  input  uj  is  given  by  CBuj^.  The  method,  proposed  by 
[8],  consists  in  finding  an  admissible  control  u  such  that  CBu  optimally  approximates  the  desired 
vaJue  CBurf.  The  optimization  criterion  is  a  least-squares  criterion,  so  that  the  objective  is  to  find 
u  €  17  such  that 

e(u,Ud)=||CBu-CBud||2  (4.10) 

is  minimized.  In  the  absence  of  constraints,  the  solution  can  be  obtained  by  setting 

^  II  CBu- CBui  IP  =  ^(CBu-CBuj)^(CBu-CBurf) 

=  2(CB)^(CBu-CBui) 

=  0.  (4.11) 

In  the  case  of  an  invertible  CB  matrix,  (4.11)  yields  u  =  uj.  With  constraints,  the  solution 
can  be  obtained  in  a  similar  way.  The  algorithm  considers  zero,  one,  two,  and  three  constraints, 
respectively,  until  a  solution  can  be  determined. 

Algorithm 

1.  If  uj  e  17,  let  u  =  U(i  and  stop.  Else,  proceed. 

2.  For  ui  =  til, mini  fiitd  u  such  that  e(u,U(i)  given  in  (4.10)  is  minimized  without  constraints 
on  U2  and  t/3.  To  that  effect,  define  U2,3  to  be  the  vector  composed  of  the  second  and  third 
elements  of  u,  (CB)i  to  be  the  vector  composed  of  the  first  column  of  CB,  and  (CB)2,3  to 
be  the  3  X  2  matrix  composed  of  the  second  and  third  columns  of  CB.  Then,  the  optimum  u 
is  given  by 

U2.3  =  [(CB)j3(CB)2.3]-'(CB)i;3[(CB)uj  -  (CB)iUi, (4.12) 
Determine  the  value  of  e(u,  u^;). 

3.  Repeat  step  2  with  the  appropriate  change  of  indices  to  determine  the  optimal  u’s  constrained 

to  til  '^il.mori  and  similarly  for  ti2  —  ti2,min»  ^2  —  ii2,max>  ii3  —  and  ti3  =  ti3,ntai' 

Determine  the  value  of  e(u,  u^)  in  every  case. 

4.  For  til  =  wi,m(n»  W2  =  W2,min,  find  u  such  that  e(u,U(i)  is  minimized  without  constraints  on 
ti3,  that  is 

ti3  =  [(CB)J(CB)3]-'(CB)J[(CB)u^  -  (CB)iUi,„,.„  -  (CB)2ti2.,n.„].  (4.13) 

Determine  the  value  of  e(u,  uj). 

5.  Repeat  step  4  with  the  appropriate  change  of  indices  to  determine  the  optimal  u’s  constrained 
to  til  =  Ui^min,  «2  =  U2,max  and  similarly  for  every  pair  among  tii  =  tii,mtn,  =  wi.max, 
'^2  =  1^2, min,  U2  =  U2,max,  W3  =  Uz^rnin,  and  U3  =  u^^^ax-  Determine  the  value  of  e(u,  uj)  in 
every  case. 

6.  Determine  the  value  of  e(u,  u^i)  for  tii  =  tii,„i,-n,  U2  =  U2,min,  ^3  =  U3,min,  and  similarly  for 
every  triple  among  tii  —  ‘^i,min,  —  '^i,max,  '^2  —  '^2, min,  1^2  —  '^2, max,  'ti3  —  f^3,min,  and 

7.  Collect  ail  the  candidate  u’s,  eliminate  those  that  do  not  belong  to  U ,  and  select  the  one  that 
yields  the  smallest  e(u,  uj). 

The  algorithm  can  be  simplified  if  a  solution  is  determined  within  the  first  steps.  For  example, 
step  4  is  not  necessary  if  an  admissible  u  is  determined  in  step  2.  The  same  is  true  in  steps  5  and  6 
if  a  solution  is  obtained  in  the  previous  steps. 
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4.6  Simulation  Results 


Figure  4.4  shows  the  results  of  applying  Methods  #1  &  #2.  Figure  4.5  shows  the  results  using 
Methods  #3  &  ^4.  Comparing  the  results  of  these  figures  with  the  results  without  command 
limiting  shown  in  Figure  4.2,  one  finds  that  responses  are  significantly  improved  by  each  of  the  four 
methods.  All  methods  perform  comparably  for  the  pitch  rate  responses  shown  in  plot  (a)  of  each 
figure.  The  ringing  is  eliminated,  and  the  coupling  is  significantly  reduced.  The  coupling  continues 
to  decay  with  time  in  the  roll  rate  responses  shown  in  plot  (b)  of  each  figure,  unlike  the  response 
shown  in  Figure  4.2  (b).  Furthermore,  performance  is  comparable  with  regard  to  coupling  for  all 
methods,  although  some  transient  oscillations  appear  with  Method  #  4.  The  similar  performance 
of  the  four  methods  is  surprising  since  their  principle  and  their  relative  complexity  are  so  different 
from  each  other.  A  possible  explanation  is  that  all  four  methods  manage  to  restore  the  directionality 
of  the  commands  to  a  reasonable  level,  while  simple  saturation  of  the  controls  fails  to  do  so. 

The  commands  were  limited  about  20%  of  the  time  in  the  pitch  rate  command  runs.  For 
Method  #2,  the  “no  p”  case  was  encountered  2%  of  the  time,  i.e.,  10%  of  the  time  the  command 
limiting  routine  was  called.  For  Method  #3,  the  “no  p”  case  was  encountered  about  0.4%  of  the 
time. 


4.7  Summary 

The  adaptive  control  algorithm  had  shown  significantly  degraded  performance  for  larger  step  size 
maneuvers,  partially  due  to  saturation  occurring  in  the  actuators.  The  control  law  would  calcu¬ 
late  command  values  that  were  not  attainable  by  the  actuators.  The  problem  was  alleviated  by 
developing  four  different  methods  of  command  limiting.  Simulation  results  do  not  show  one  par¬ 
ticular  method  to  be  superior  to  the  others,  although  the  complexity  is  quite  different  for  the  four 
methods.  All  four  methods  do,  however,  offer  significant  improvement  to  the  performance  of  the 
aircraft.  To  simplify  future  analyses,  the  least-squares  approximation  to  commanded  accelerations 
will  be  used. 
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Rol  Rato  Response  to  Pitch  Rate  Command:  Model.  Methods  1  &  2 


(a)  (b) 


Right  Stabilalor  Position  for  Pitch  Rate  Command;  Methods  1  &  2  Right  Aileron  Position  for  Pitch  Rate  Command:  Methods  1  &  2 


—  model  response 

- Method  95^1 
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Figure  4.4:  Responses  to  a  Large  Pitch  Rate  Command  Using  Methods  #1  &  #2 
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Time  (sec)  Time  (sec) 

(c)  (d) 
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Figure  4.5:  Responses  to  a  Large  Pitch  Rate  Command  Using  Methods  &  #4 
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Chapter  5 

Steady- State  Tracking 


Integral  compensation  is  used  in  nonadaptive  linear  control  systems  to  ensure  the  tracking  of  con¬ 
stant  reference  inputs  and  the  rejection  of  constant  disturbances.  Since  the  responses  of  the  aircraft 
using  the  adaptive  flight  controller  have  shown  small  errors  in  steady-state  tracking,  particularly 
in  times  of  failure,  the  application  of  integral  compensation  to  the  control  law  was  considered. 
Integral  compensation  can  be  incorporated  in  the  adaptive  control  algorithms  discussed  so  far.  The 
improvements  that  can  be  obtained  with  this  addition  will  be  shown. 


5.1  Benefits  of  Integral  Control 

Integral  compensation  is  reviewed  in  general  terms  in  this  section.  A  single-input,  single-output 
plant  P{s)  is  considered  with  input  u,  output  y,  and  input  disturbance  d,  so  that 

y  =  P{s)[u -{■  d].  (b.l) 

The  term  d  represents  a  constant  disturbance,  i.e.,  d(t)  -  d”.  It  was  shown  in  Section  2.2  that 
the  trim  uncertainty  could  be  viewed  as  a  constant  disturbance,  and  a  major  objective  of  the 
reconflgurable  control  system  is  to  compensate  for  trim  variations.  A  standard  compensator,  based 
on  the  tracking  error  tm  —  y,  is  given  by 


u  =  C(s)[rM  -  y], 

where  tm  is  the  reference  input.  The  resulting  output  is 

y  =  H,{s)[rM]^  H2{s)[d], 


(5.2) 


(5.3) 


where 


^i(s) 


P{s)C{s) 

i  +  p{s)c{sy 


this) 


Pjs) 

i  +  p{s)c{sy 


(5.4) 


Integral  compensation  consists  in  having  a  pole  of  C(s)  at  s  =  0.  In  that  case,  Ifi(O)  =  1,  which 
means  that  constant  reference  inputs  are  perfectly  tracked,  and  H^iO)  =  0,  which  means  that 
constant  disturbances  are  perfectly  rejected.  Attractive  features  of  integral  compensation  are  that 
both  of  the  objectives  are  achieved  together,  and  that  they  are  achieved  despite  uncertainties  in 
the  plant  parameters,  as  long  as  the  closed-loop  system  remains  stable. 


35 


Figure  5.1:  Model  Reference  Control  Loop  with  Integral  Compensation 

5.2  Integral  Control  Law 

In  [18],  Bodson  and  Pohlchuck  proposed  to  add  integral  compensation  to  the  control  law.  The 
original  control  law  of  (2.12)  is  shown  in  Figure  2.7.  The  modified  control  law  is 


u  =  Co(rAr  -  y  +  ^-  /  (tm  -  y)dt)  +  Gqx  +  vq  (5.5) 

Jo 

and  is  shown  in  Figure  5.1.  Assuming  zero  initial  conditions  and  the  aircraft  model  y=  CAx + 
CBu  +  Cd,  one  has,  in  the  Laplace  domain 

(s^  +  2ks  +  k^)'y{s)  =  CAsx(s)  +  CBsu(s)  +  Csd(s)  +  2ksy{s)  +  k^y{s) 

=  (CA  +  CBGo)sx(s)  +  CBCo5rM(s)  -  CBCo5y(s) 
+A:CBCorA/(s)  -  A:CBCoy(5)  +  CB5Vo(5)  +  Csd(s) 
+2ksy{s)  +  k'^y{s).  (5.6) 

Using  the  nominal  values  of  the  parameters  yields 

In  other  words,  using  the  same  parameter  matrices,  the  closed-loop  response  remains  the  first-order 
reference  model  transfer  function  of  (2.9). 


5.3  Simulation  Results 

Figures  5.2  -  5.5  compare  results  with  and  without  integral  compensation.  The  least-squares 
command  limiting  method  is  applied.  Responses  with  a  stabilator  locked  at  305  seconds  during 
a  batch  least-squares  run  appear  in  Figure  5.2.  The  addition  of  integral  compensation  improves 
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Right  SUbilaUx  Position  (deg)  q  (deg/s) 


Pitch  Rate  Response  to  Pitch  Rate  Command:  Model.  BLS,  &  IBLS 


Roll  Rate  Response  to  Pitch  Rata  Command:  Modal.  BLS.  &  IBLS 


Time  (sec)  Time  (sec) 

(c)  (d) 

—  model  response 

- no  integral  compensation 

—  •  —  integral  compensation 

Figure  5.2:  BLS  Results  for  the  Locked  Stabilator  Aircraft  with  and  without  Integral  Compensation 


37 


steady-state  tracking  and  prevents  drift  in  this  case.  Interestingly,  improvements  in  cross-coupling 
effects  are  also  observed  in  plot  (b).  Thus,  for  the  batch  least-squares  case,  integral  compensation 
generates  near-model  responses  in  both  the  pitch  and  roll  rate. 

Results  are  shown  in  Figure  5,3  for  the  stabilized  recursive  least-squares  control  algorithm  ap¬ 
plied  to  an  aircraft  in  which  the  left  stabilator  becomes  locked  at  310  seconds.  Note  that  the 
tracking  errors  are  slightly  smaller  with  the  recursive  least-squares  algorithm  than  with  the  batch 
least-squares  algorithm.  Furthermore,  integral  compensation  causes  the  response  of  the  aircraft  to 
more  quickly  match  the  model  response  once  a  failure  occurs.  Observing  the  right  stabilator  and 
aileron  positions  in  plots  (c)  and  (d),  respectively,  shows  that  after  the  failure  occurs,  larger  com¬ 
mands  to  both  actuators  result.  The  trim  to  both  actuators  also  changes,  and  integral  compensation 
prevents  divergence  of  the  aileron  position. 

Figure  5.4  considers  the  recursive  least-squares  control  algorithm  applied  to  the  aircraft  with 
a  missing  left  stabilator  occurring  at  310  seconds.  Again,  integral  compensation  causes  the  model 
response  to  be  matched  more  quickly.  Similar  to  the  previous  case,  the  cross-couplings  in  plot  (b) 
are  lower  with  integral  compensation.  Note,  however,  a  period  of  rapid  oscillation  occurs  with 
integral  compensation. 

Figure  5.5  considers  the  locked  stabilator  failure  for  a  step  size  of  ±10®/s.  Integral  compensation 
does  not  adversely  affect  the  aircraft  during  this  maneuver.  However,  it  does  not  seem  to  offer  great 
improvement  either.  In  fact,  integral  compensation  generates  a  higher  transient  cross-coupling 
in  plot  (b).  Overall,  integral  compensation  appears  to  significantly  improve  batch  least-squares 
responses,  but  only  offers  slight  improvement  in  the  recursive  least-squares  cases. 

5.4  Summary 

Improvements  to  steady-state  tracking  can  be  achieved  by  adding  integral  compensation  to  the 
control  law.  The  most  significant  improvement  results  using  the  batch  least-squares  identification 
algorithm.  This  is  compatible  with  conventional  knowledge  regarding  the  use  of  integral  compensa¬ 
tion  in  nonadaptive  systems.  Integral  compensation  does  not  offer  much  improvement  when  using 
the  recursive  least-squares  identification  algorithm,  i.e.,  in  the  adaptive  context  under  considera¬ 
tion  for  this  application.  The  recursive  algorithm  is  able  to  fine-tune  the  controller  matrices  and, 
as  a  result,  provide  some  of  the  benefits  normally  attributed  to  integral  compensation.  Addition¬ 
ally,  integral  compensation  using  the  recursive  least-squares  algorithm  often  causes  the  parameters 
to  change  more  abruptly,  thus  resulting  in  spiked  or  oscillatory  behavior  in  the  aircraft  output. 
Overall,  the  benefits  of  integral  compensation  are  limited  and  probably  not  worth  the  additional 
complexity. 
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RJghI  Slabilalof  Portion  (deg)  q  (deg/s) 


Figure  5.3:  RLS  Results  for  the  Unfailed  then  Locked  Stabilator  Aircraft  with  and  without  Integral 
Compensation 
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—  model  response 

- no  integral  compensation 
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Figure  5.4:  RLS  Results  for  the  Unfailed  then  Missing  Stabilator  Aircraft  with  and  without  Integral 
Compensation 
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- no  integral  compensation 
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Figure  5.5:  RLS  Results  for  the  Unfailed  then  Locked  Stabilator  Aircraft  for  a  Large  Command 
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Chapter  6 

Noise  Analysis 


6.1  Noise  Models 

Measurement  devices  generate  signals  that  contain  noise.  The  addition  of  noise  can  complicate 
the  control  problem,  especially  if  adaptation  is  required.  Previous  work  did  not  include  the  effect 
of  noise  because  noise  models  were  not  part  of  the  original  simulation  code.  However,  concern 
about  the  possible  impact  of  the  noise  in  a  practical  environment  needed  to  be  addressed.  For 
this  reason,  a  few  experiments  were  performed  with  noise.  Noise  was  modeled  as  Gaussian  white 
noise  with  standard  deviations  determined  from  the  literature  [19]  as  well  as  flight  test  data,  and 
was  added  to  the  measurements  used  by  the  adaptive  algorithms  and  the  control  law.  Table  6.1 
lists  the  standard  deviations  used  in  the  simulations  for  various  signals.  In  the  code,  a  Gaussian 
noise  source  is  created  from  a  uniform  random  number  generator  [20],  [21]  converted  to  a  Gaussian 
distribution  [22].  The  specific  routines  for  generating  the  Gaussian  random  numbers  are  described 
in  Appendix  D. 

6.2  Simulation  Results 

Figures  6.1  and  6.2  illustrate  the  performance  of  the  adaptive  algorithm  with  and  without  the  noise 
models  applied  for  step  sizes  of  ±l°/s  and  ±10°/s,  respectively.  The  algorithm  uses  the  least- 


Table  6.1:  Noise  Standard  Deviations 


Parameter 

Noise  Standard  Deviation 

a 

0.1° 

/? 

0.1° 

9 

0.1°/s 

V 

0.1°/s 

r 

0.1°/s 

9 

1.5°/s2 

V 

3.0°/s2 

r 

1.5°/s''^ 
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squares  approximation  of  commanded  accelerations  as  the  method  of  command  limiting.  Results 
illustrate  that  the  algorithm  handles  the  noise  well.  Results  appear  to  be  worse  in  the  ±l°/s  case 
because  the  values  of  the  noise  are  similar  in  magnitude  to  the  step  sizes  commanded.  The  noise 
does  not  appear  to  have  much  of  an  effect  on  the  ±10°/s  run.  These  runs  are  representative  of  the 
behavior  of  the  adaptive  control  law  observed  for  a  number  of  simulations  in  which  noise  was  added. 
No  simulations  with  noise  generated  unstable  behavior  in  the  aircraft.  Overall,  the  algorithm  was 
found  to  be  capable  of  handling  the  noise  well. 

As  previously  mentioned,  the  method  of  command  limiting  used  in  Figures  6.1  and  6.2  was  the 
least-squares  approximation  of  commanded  accelerations.  Simulations  were  also  run  with  the  other 
command  limiting  methods.  Interestingly,  the  noise  affected  the  amount  of  command  limiting.  For 
example,  in  Figures  4.4  and  4.5  commands  were  limited  about  20%  of  the  time.  When  noise  was 
added  to  the  same  simulation,  the  commands  were  limited  around  25%  of  the  time.  In  Figures  4.4 
and  4.5  the  “no  p"  case  was  encountered  about  2%  of  the  time.  With  the  addition  of  noise, 
“no  p”  cases  occurred  about  5%  of  the  time.  However,  the  increase  in  command  limiting  did  not 
significantly  reduce  the  performance  of  the  adaptive  controller. 


6.3  Summary 

The  brief  analysis  of  this  chapter  showed  that  noise  did  not  affect  the  adaptive  algorithm  signif¬ 
icantly.  A  simple  Gaussian  white  model  of  noise  added  to  the  sensors  was  used.  In  practice,  the 
perturbation  observed  in  the  measured  signals  consists  primarily  of  quantization  and  unmodeled 
effects  such  as  flexible  modes.  This  noise  would  be  more  complicated  to  model  and  would  be  the 
focus  of  future  research.  The  above  results,  however,  demonstrate  that  the  algorithm  is  not  unduly 
sensitive  to  realistic  amounts  of  measurement  noise. 
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Right  Stabilatof  Position  (deg)  q  (deg/s) 


(a) 


(b) 


flight  Stabilalor  Positbn  for  Pitch  Rate  ComrnarKf:  w/o  &  w/  Noisa  Right  Aileron  Position  for  Pitch  Rata  Command:  w/o  &  w/  Noise 


—  model  response 

- without  noise 

—  •  —  with  noise 


Figure  6.2:  Results  with  and  without  Noise  for  ±10®/s  Locked  Stabilator  Failure 


Chapter  7 

Evaluation  for  Different  Flight 
Conditions  and  Failures 


This  chapter  examines  the  performance  of  the  adaptive  controller  for  other  flight  conditions,  in¬ 
cluding  two  other  subsonic  and  one  supersonic  condition.  Additional  failures,  a  locked  aileron  and 
a  combined  locked  stabilator  and  aileron,  are  also  examined. 

7.1  Flight  Condition:  9,800  feet,  Mach  0.9 

Figure  7.1  illustrates  the  performance  of  the  aircraft  for  flight  condition  9,800  feet,  Mach  0,9,  for 
a  locked  stabilator.  The  response  of  the  aircraft  is  comparable  to  the  previous  responses,  although 
the  tracking  is  somewhat  degraded.  A  similar  degradation  is  observed  for  the  39,800  feet,  Mach  1.4 
condition,  and  its  discussion  is  postponed  until  later. 

7.2  Flight  Condition:  39,800  feet,  Mach  0.6 

The  39,800  feet,  Mach  0.6  flight  condition  requires  a  slight  modification  to  the  testing  procedure. 
As  mentioned  in  Section  2.5  and  Appendix  A,  the  proportional-integral  controller  is  used  to  fly 
the  aircraft  for  the  first  305  seconds  of  flight  in  order  for  transients  to  die  out.  Certain  aircraft 
parameters,  specifically  the  altitude  and  airspeed,  must  be  set  during  this  time.  For  flight  conditions 
well  inside  the  operational  envelope,  setting  the  altitude  and  airspeed  to  the  values  specified  by  the 
flight  condition  does  not  pose  a  problem  using  the  proportional-integral  controller.  However,  the 
flight  condition  of  39,800  feet,  Mach  0.6  is  close  to  the  boundary  of  the  operational  envelope  [7]. 
If  the  altitude  and  airspeed  are  set  to  the  flight  condition  values  in  this  case,  the  proportional- 
integral  controller  does  not  maintain  flight  for  305  seconds.  Hence,  by  305  seconds  the  condition 
of  the  aircraft  is  not  39,800  feet,  Mach  0.6.  The  slight  modification  incorporated  in  the  code  was 
to  adjust  the  set  point  of  the  airspeed.  The  airspeed  is  set  to  678  ft/s  (Mach  0.7)  for  the  first 
200  seconds  of  flight  and  to  631  ft/s  (Mach  0.65)  afterwards.  It  may  be  noted  that  631  ft/s  does 
not  equal  Mach  0.6,  but  in  this  case  the  set  point  generates  an  actual  aircraft  speed  of  Mach  0.6. 

Figure  7.2  shows  the  response  of  the  aircraft  to  a  locked  stabilator  failure.  The  performance 
of  the  aircraft  is  still  reasonable  but  is  degraded  under  this  flight  condition.  Cross-couplings  in 
plot  (b)  are  higher  in  response  to  the  pitch  rate  command,  and  the  pitch  rate  response  in  plot  (a) 
is  delayed.  The  degradation  is  related  to  the  extensive  rate  saturation  shown  by  the  linear  portions 
of  the  right  stabilator  and  aileron  positions  shown  in  plots  (c)  and  (d),  respectively.  At  this  flight 
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Righl  SlabiUtof  Position  (d«g)  q  (dag/s) 


10, 


Pitch  Rata  Rasponsa  to  Pitch  Rala  Command:  9800  feet  Mach  0.9 


6i 


RoO  Rale  Rasponsa  to  Pitch  Rala  Command:  9800  feat  Mach  0.9 


—  model  response 
- aircraft  response 

Figure  7.1:  Responses  for  a  Flight  Condition  of  9,800  feet  and  Mach  0.9 
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Right  SUbllator  Position  (deg)  q  (deg/s) 


Pitch  Rate  Response  to  Pilch  Rate  Command:  39800  (set  Mach  0,6 


Roll  Rate  Response  to  Pitch  Rata  Command;  39800  (set  Mach  0.6 


(a) 


(b) 


Right  Stabitatof  Position  lor  Pitch  Rate  Command;  39800  leet  Mach  0.6  Right  Aileron  Position  (or  Pitch  Rale  Command:  39800  feel  Mach  0.6 


(c) 


(d) 


—  model  response 
- aircraft  response 


Figure  7.2:  Responses  for  a  Flight  Condition  of  39,800  feet  and  Mach  0.6 


condition,  the  eflFectiveness  of  the  actuators  is  much  reduced  because  of  the  higher  altitude  and 
lower  air  density.  This  fact  can  be  observed  in  the  state-space  matrices  found  in  Appendix  A. 
The  (2, 1)  element  of  the  B  matrix  in  (2.4)  indicates  the  effectiveness  of  the  stabilator  to  generate 
a  pitching  moment.  The  greater  the  magnitude  of  B(2, 1),  the  more  effective  the  stabilator  is. 
For  the  first  flight  condition,  B(2,l)  =  -9.98.  For  the  flight  condition  of  39,800  feet,  Mach  0.6, 
B(2, 1)  =  -2.88.  The  performance  is  degraded  simply  because  a  ±107s  step  is  a  much  more 
demanding  maneuver  for  this  flight  condition  than  the  other  conditions.  Despite  this  situation,  the 
control  system  performs  well. 


7.3  Flight  Condition:  39,800  feet,  Mach  1.4 


Figure  7.3  shows  the  performance  of  the  aircraft  for  the  supersonic  flight  condition  of  39,800  feet, 
Mach  1.4  for  a  step  size  of  ±l®/s  with  and  without  integral  compensation  for  a  locked  stabilator 
failure.  A  first  observation  for  this  flight  condition  is  that  pitch  rates  are  limited,  although  not 
because  of  low  pitching  authority,  but  because  of  low  aileron  effectiveness.  A  step  size  of  ±l°/s  can 
maintain  low  coupling,  as  shown  in  plot  (b).  For  step  sizes  greater  than  ±l®/s,  the  ailerons  saturate 
in  position.  Even  though  the  effectiveness  for  pitching  is  much  higher  at  this  flight  condition  than 
at  the  nominal  condition,  the  rolling  effectiveness  is  considerably  lower,  making  it  difficult  to  pitch 
the  aircraft  without  rolling. 

Another  observation  is  that  tracking  of  the  pitch  rate  command  is  not  good  without  integral 
compensation,  although  none  of  the  actuators  saturate  in  any  way.  Simulations  of  the  unfailed  air¬ 
craft  using  a  BLS  identification  show  a  similar  behavior  with  a  steady-state  value  of  only  ±0.7®/s 
for  a  commanded  step  of  il°/s.  The  problem  was  found  to  be  due  to  the  actuator  dynamics,  and 
could  be  resolved  in  different  ways,  one  of  them  being  integral  compensation.  If  integral  compen¬ 
sation  is  applied  to  the  unfailed  aircraft,  the  steady-state  value  of  ±l°/s  is  achieved.  Figure  7.3 
shows  that  integral  compensation  improves  the  steady-state  tracking  after  the  failure. 

Unlike  the  previous  flight  conditions  in  which  the  differential  stabilator  could  be  enabled,  i.e.,  by 
letting  Suf.  =  enabling  the  differential  stabilator  in  the  supersonic  case  yields  very  poor  control 
of  the  aircraft.  The  poor  performance  is  indicated  by  the  state-space  matrices.  The  B  matrix  was 
obtained  for  the  locked  stabilator  aircraft  at  the  supersonic  flight  condition  with  6dc  0.  CB  is 
given  by 


CB  = 


-6.65 

-19.59 

-1.47 


5.58 

19.36 

1.22 


-0.20 

-1.93 

-3.45 


(7.1) 


The  first  and  second  columns  of  CB  are  the  coefficients  of  6hc  and  respectively.  The  approx¬ 
imate  equal  ratio  (~  3)  between  the  first  two  rows  in  the  first  and  second  columns  indicates  the 
moments  produced  for  pitch  and  roll  commands  are  nearly  identical  in  direction.  The  adaptive 
algorithm  cannot  compensate  for  the  failure  under  this  condition.  However,  with  again  setting 
6dc  —  0,  CB  becomes 


CB  = 


-6.90 

0.03 

-0.10 

-20.04 

2.97 

-1.49 

-1.69 

0.02 

-3.47 

(7.2) 


Although  this  matrix  is  stiU  poorly  conditioned,  the  simulation  results  were  much  better  with  this 
design. 
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Right  Stabilatof  Position  (deg)  q  (deg/s) 


Pitch  Rata  Response  to  Pitch  Rate  Command;  39800  (eet  Mach  1 .4  RoB  Rate  Response  to  Pitch  Rate  Command;  39800  feet  Mach  1.4 


Ttme  (sec)  Time  (sec) 

(c)  (cl) 

—  model  response 

- response  without  integral  compensation 

—  •  —  response  with  integral  compensation 

Figure  7.3:  Responses  for  a  Flight  Condition  of  39,800  feet  and  Mach  1.4 
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7.4  Locked  Aileron  Failure 


The  adaptive  algorithm  is  effective  for  failures  other  than  the  locked  stabilator  and  missing  stabi- 
lator.  A  locked  aileron  failure  is  considered  for  illustration.  As  opposed  to  the  stabilator  failures 
which  induce  rolling  moments,  the  aileron  failure  does  not  induce  significant  pitching  moments, 
so  that  this  failure  is  less  demanding.  Figures  7.4  and  7.5  illustrate  the  effectiveness  of  the  com¬ 
mand  limiting  techniques  to  roll  rate  commands  with  step  sizes  of  ±30®/$  and  ±50°/s,  respectively. 
For  the  responses  in  Figure  7.4  for  ±30°/s,  very  little  difference  is  found  between  the  techniques. 
Couplings  in  plots  (a)  and  (c)  remain  low,  and  the  response  nearly  approaches  the  model  in  plots  (b) 
and  (d).  Command  limiting  was  applied  about  20%  of  the  time,  and  the  “no  p”  case  occurred  less 
than  1%  of  the  time.  The  responses  in  Figure  7.5  for  ±50®/s  do  not  provide  the  same  conclu¬ 
sions.  The  command  limiting  methods  perform  well,  except  the  least-squares  approximation  to 
commanded  accelerations,  i.e..  Method  #  4.  The  least-squares  method  has  an  extremely  large 
coupling  that  develops  for  a  period  of  time,  although  the  response  is  recovered  as  shown  in  plots  (c) 
and  (d).  For  the  step  size  of  ±50®/s,  command  limiting  occurred  about  35%  of  the  time,  except  in 
the  least-squares  case  in  which  command  limiting  occurred  about  60%  of  the  time.  The  significant 
amount  of  command  limiting  can  be  related  to  the  transient  response  observed  in  Figure  7.5.  For 
the  relaxation  of  control  requirements  and  the  scaling  of  reference  inputs  methods,  the  “no  p”  case 
occurred  less  than  3%  of  the  time. 

7.5  Combined  Stabilator/ Aileron  Failure 

A  combined  stabilator  and  aileron  failure  was  also  applied  to  the  aircraft  model.  In  this  case,  both 
the  left  stabilator  and  left  aileron  were  locked  at  310  seconds.  Figure  7.6  shows  the  performance  of 
the  aircraft  with  the  combined  failure  to  a  pitch  rate  command  of  ±10®/s.  The  control  algorithm  is 
able  to  achieve  a  satisfactory  response.  The  deviations  of  both  the  unfailed  stabilator  and  aileron, 
shown  in  plots  (c)  and  (d),  respectively,  increase  significantly  due  to  this  failure.  The  algorithm  is 
also  capable  of  handling  a  roll  rate  command  of  ±30°/s  as  shown  in  Figure  7.7. 

7.6  Summary 

The  adaptive  algorithm  was  shown  to  perform  well  for  other  flight  conditions  and  failure  types. 
Degradation  in  performance  at  various  flight  conditions  was  due  to  a  variety  of  factors  including 
loss  of  actuator  effectiveness.  At  one  flight  condition,  actuator  dynamics  had  a  significant  impact. 
Further  research  is  needed  to  address  this  problem.  The  adaptive  controller  has  been  shown  to  be 
capable  of  controlling  the  aircraft  during  a  locked  aileron  and  a  combined  locked  stabilator  and 
locked  aileron  failure.  The  adaptive  algorithm  generally  performed  well  for  large  step  sizes  in  roll 
rate,  with  roll  rate  most  significantly  affected  by  an  aileron  failure.  Some  interesting  behavior  that 
resulted  was  that  the  least-squares  approximation  to  commanded  accelerations  seemed  to  perform 
worse  for  a  very  large  roll  rate  command  with  a  locked  aileron. 
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Pitch  Rata  Response  to  Roll  Rata  Command:  Modal.  Methods  1  &  2 
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Roll  Rale  Response  to  Roll  Rale  Command;  Model,  Methods  1  &  2 
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Figure  7.4:  Roll  Rate  Command  Responses  for  ±30°/s  and  All  Command  Limiting  Methods  with 
a  Locked  Aileron  Failure 
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Figure  7.5:  Roll  Rate  Command  Responses  for  ±50°/s  and  All  Command  Limiting  Methods  with 
a  Locked  Aileron  Failure 


15| 


Pitch  Ralo  Raaponse  to  Pitch  Rate  Command:  Combined  Faiturea 


RoD  Rate  Response  to  Pitch  Rate  Command;  Combined  Failures 


—  model  response 
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Figure  7.6:  Pitch  Rate  Command  Responses  for  ±10®/s  with  a  Combined  Locked  Stabilator/ Aileron 
Failure 
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Right  Stabilalor  Po&ibon  (dag)  q  (dag/s) 


Pitch  Rate  Response  to  Roll  Rale  Command;  Combined  Failures 


Ron  Rate  Response  to  Roll  Rale  Command:  Combined  Failures 


—  model  response 
- aircraft  response 


Figure  7,7:  Roll  Rate  Command  Responses  for  ±30® /s  with  a  Combined  Locked  Stabilator/ Aileron 
Failure 
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Chapter  8 

Miscellaneous  Observations 


A  few  observations  made  in  the  course  of  the  research  were  besides  the  focus  of  the  research.  This 
chapter  addresses  these  issues. 


8.1  Slow  Divergence 

The  adaptive  controller  does  not  regulate  pitch,  roll,  and  yaw  angular  positions  in  the  control  law. 
Since  these  values  are  the  integrals  of  the  rates,  divergence  may  easily  result.  In  practice,  the  pilot 
would  regulate  angular  position  given  control  of  the  angular  rates.  In  the  simulation,  an  autopilot 
may  be  applied  to  return  the  aircraft  to  nominal  values  of  pitch,  roll,  and  yaw  angles  and  prevent 
the  divergence.  The  autopilot  matches  the  model  (2.8),  (2.10),  (2.11)  and  is  given  by 

9  =  -2.5(7  +  2.5(7c, 

i>  =  -2.5p  +  2.5pc,  (8.1) 

r  =  — 2.5r  +  2.5rc. 

Using 

0  =  9, 

4>  =  p,  (8.2) 

i)  =  r, 

where  0,  (f>,  and  V’  are  the  angular  positions  of  pitch,  roll,  and  yaw,  respectively,  and  using  the 
commands 


qc  ~  di^c  ^)> 

Pc  =  (8.3) 

re  =  oiA  -  V’), 

the  autopilot  will  track  the  desired  angles  with  the  transfer  function 


2.5g 

"h  2.55  2.5;^ 


(8.4) 


The  constant  g  is  selected  to  be  1.6  so  that  the  closed-loop  poles  are  located  at  —1.25  ±  71.56. 
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PHch  Rale  Response  to  PHch  Rale  Command:  BLS  w/o  &  w/  Autopilot  Correction 


(a) 


Rotl  Rato  Response  to  Pitch  Rate  Command:  BLS  w/o  &  w/  Aulopilol  Correction 


. ; . i'"' 

1  ' 

.»  •  1  ' 

. .'  'i'' . ;  'r 

' 

. / . i- .  1 

. »  V.’ . t 

! 

. . ; . 

1  'n . i . . 

'  t  A  f 

. f . . ' 

''  »  •  '  1 
\  1  •  *  ' 

*,  •  :  '  1 

1  '  ^  ' 

C _ 1 _ _ : _ 1 _ ^ _ 

-  r  1  i 

:  \ 

- 1 _ 1 _ 

- - 1 - 

t 

_ 1 

i - \ - 1 _ I 

325  330  335  340  345 

Time  (sec) 
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—  model  response 

- no  autopilot  correction 

—  •  —  with  autopilot  correction 


Figure  8.1:  Demanding  Flight  Maneuver  Rate  Responses  with  and  without  Autopilot  Correction 

Figure  8.1  compares  the  unfailed  batch  least-squares  responses  with  and  without  applying  an 
autopilot  correction  in  between  stepping  segments  of  size  ±10°/s  without  control  allocation.  Drift 
in  the  roll  rate  response  to  a  pitch  rate  command  shown  in  plot  (b)  is  observed  without  the 
autopilot.  Overall,  the  responses  show  that  slow  divergence  due  to  the  uncontrolled  states  can 
easily  be  compensated  for. 

8.2  Glitches  in  Responses 

Some  simulations  show  behavior  that  is  unexpected.  Figure  8.2  shows  a  simulation  for  the  unfailed 
aircraft  where  the  length  of  the  steps  is  extended  to  5  seconds,  as  opposed  to  1.25  seconds.  As 
the  pitch  rate  response  approaches  the  final  value  of  —  l°/s,  a  kink  is  observed  in  the  response  of 
plot  (a),  which  makes  the  response  look  like  a  double  step  response.  Such  glitches  often  appear 
in  the  simulations  and  are  due  to  the  discrete  nature  of  look-up  tables  used  in  the  simulations. 
Specifically,  Figure  8.2  shows  the  index  of  the  look-up  tables  based  on  the  angle  of  attack.  It  is 
superimposed  on  the  pitch  rate,  roll  rate,  and  angle  of  attack  plots  and  is  scaled  for  easy  comparison. 
The  index  is  IDXAAl  and  is  shown  as  a  dashed  line  in  the  figure.  IDXAAl  changes  for  every  4° 
change  in  the  angle  of  attack.  Viewing  plots  (a)  and  (c)  of  the  figure  shows  that  the  glitch  occurs 
at  a  point  where  the  angle  of  attack  index  changes.  The  glitch  is  not  a  problem  of  the  controller 
but  an  imperfection  in  the  interpolation  of  the  look-up  tables. 

8.3  Nonlinear  Effects  of  the  Angle  of  Attack 

Extending  the  length  of  time  for  the  pitch  rate  command  produces  a  problem  for  the  control  law, 
because  of  large  angles  of  attack.  The  standard  length  of  time  for  a  step  is  1.25  seconds.  Extending 
the  step  time  length  to  2.5  seconds  causes  a  larger  excursion  in  angle  of  attack.  A  ±10°/s  step  in 
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Pitch  Rata  Responaa  to  Pitch  Rata  Command:  Step  Length  o(  S  aec 


Roll  Rate  Reaporwa  to  Pitch  Rate  Command:  Step  Length  of  5  $ec 


(a)  (b) 


Angle  of  Attack  and  tOXAAl  for  a  Pitch  Rata  Command:  Step  Length  of  5  sec 
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—  aircraft  response 
- scaled  IDXAAl  index 


Figure  8,2:  Unfailed  Aircraft  Responses  to  a  Pitch  Rate  Command  for  Step  Time  Length  of  5  sec 


pitch  causes  the  angle  of  attack  to  change  from  4°  to  about  20°  in  the  2.5  second  interval.  Figure  8.3 
illustrates  that  spikes  occur  in  the  coupling  response  in  plot  (b)  when  the  positive  steady-state  value 
in  the  desired  axis  is  approached.  At  this  time,  the  angle  of  attack  is  approaching  20°.  The  problem 
occurring  here  is  that  the  ailerons  begin  losing  effectiveness  rapidly  as  the  angle  of  attack  increases. 
The  solution  to  this  problem  is  not  simple.  Modifications  to  the  control  law  are  required  and  would 
be  the  subject  of  future  research  in  regards  to  the  model  reference  adaptive  controller. 

8.4  Summary 

Certain  interesting  behaviors  occurred  in  the  simulations  that  were  not  directly  related  to  the 
adaptive  controller  or  were  beyond  the  scope  of  the  research.  Responses  of  the  aircraft  could  become 
poor  due  to  slow  divergence  of  the  angular  positions  of  the  aircraft,  but  this  divergence  would  be 
corrected  by  the  pilot  in  practice,  or  could  be  eliminated  by  an  autopilot.  Another  behavior  is  that 
of  glitches  appearing  in  certain  simulations  caused  by  the  aircraft  simulation  software.  The  most 
significant  problem  that  may  require  more  research  is  that  of  determining  corrective  actions  for  the 
nonlinear  effects  of  the  angle  of  attack. 
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Right  Stablialor  Position  (deg)  q  (deg/s) 


Chapter  9 

Conclusions 


The  goal  of  this  research  was  to  improve  the  robustness  of  an  adaptive  flight  control  law.  An 
algorithm  was  selected  that  was  simple,  was  capable  of  handling  the  cross-couplings  of  a  failed 
aircraft,  and  utilized  a  state- feedback  formulation.  The  basic  algorithm  combining  a  model  reference 
control  law  and  a  stabilized  recursive  least-squares  identification  procedure  was  shown  to  be  capable 
of  providing  control  reconfiguration,  but  showed  insufficient  robustness  properties. 

Two  alternate  identification  methods  were  developed  that  improved  the  performance  of  the 
control  law.  The  inverse  controller  formulation  was  a  direct  approach  that  used  the  same  matrices 
as  the  control  law,  but  required  the  calculation  of  a  matrix  inverse  to  determine  the  parameter 
matrices.  The  indirect  controller  formulation  was  similar  to  the  inverse  controller  formulation, 
with  the  added  benefit  of  removing  the  model  parameters  from  the  identification.  Instability  was 
also  eliminated  with  the  filtering  of  actuator  commands  and  normalization  of  the  algorithm.  The 
use  of  actuator  measurements  in  the  identification  slightly  improved  the  results,  but  this  option 
was  not  retained  because  additional  fault  detection  logic  would  be  required. 

A  significant  degradation  of  performance  occurred  with  larger  step  size  maneuvers,  due  to 
saturation  in  the  actuators.  To  alleviate  the  problem,  four  different  methods  of  control  allocation 
were  proposed  and  evaluated.  Simulations  did  not  show  one  particular  method  to  be  superior  to 
the  others  despite  the  differing  levels  of  complexity  of  the  algorithms.  Performance  of  the  adaptive 
algorithm  was  significantly  improved  by  all  four  methods. 

Integral  compensation  was  added  to  the  control  algorithm  to  improve  tracking  of  reference  in¬ 
puts.  However,  integral  compensation  did  not  always  offer  significant  improvement  to  the  recursive 
algorithm  because  the  recursive  algorithm  was  able  to  fine-tune  the  controller  matrices.  In  addition, 
a  recursive  algorithm  combined  with  integral  compensation  often  caused  the  parameters  to  change 
abruptly,  resulting  in  spiked  behavior.  The  benefits  of  integral  compensation  were  considered  to  be 
limited  and  possibly  not  worth  the  additional  complexity. 

A  brief  noise  analysis  was  performed.  The  results  demonstrated  that  the  algorithm  was  capable 
of  handling  Gaussian  white  noise  added  to  the  sensors.  Gaussian  white  noise  is  not  an  accurate 
model  of  sensor  noise,  but  the  results  demonstrated  that  the  algorithm  was  not  unduly  sensitive  to 
measurement  noise. 

The  adaptive  algorithm  performed  well  for  other  flight  conditions  and  failures.  Degradation  at 
some  flight  conditions  was  due  to  reduction  of  actuator  effectiveness.  In  some  other  cases,  further 
developments  of  the  algorithm  may  be  needed  to  improve  performance. 

Although  results  of  this  thesis  apply  to  a  specific  aircraft  model,  some  of  the  conclusions  should 
apply  to  other  aircrafts,  or  even  to  other  multivariable  control  problems.  First,  using  the  modified 
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indirect  identification  formulation  yields  a  much  more  stable  system  than  using  the  direct  identifi¬ 
cation  formulation.  The  indirect  identification  requires  slightly  more  calculations,  i.e.,  calculating 
Co,  Go,  and  vo,  but  the  number  of  parameters  that  are  identified  is  the  same  using  either  the  direct 
or  modified  indirect  formulation.  Second,  the  thesis  has  shown  that  rate  limiting  of  the  actuators 
can  significantly  degrade  the  performance  of  the  adaptive  control  law.  Ensuring  that  the  control 
inputs  are  modified  to  satisfy  the  rate  constraints  brings  the  performance  of  the  aircraft  to  within 
reasonable  levels.  Third,  integral  compensation  in  a  control  law  may  be  used  to  minimize  steady- 
state  tracking  errors,  but  with  less  benefits  than  might  be  expected  since  an  adaptive  algorithm 
with  bias  cancellation  is  able  to  provide  for  trim  and  steady-state  tracking.  In  addition,  spiked 
behavior  can  occur  from  the  integral  compensation. 
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Appendix  A 

Aircraft  Parameter  Matrices 


A.l  Off-Line  Identification  of  the  Parameter  Matrices 

The  indirect  parameter  matrices  A,  B,  and  d  and  the  direct  parameter  matrices  Co,  Go,  and  vo 
are  obtained  from  a  run  using  the  proportional-integral  controller.  The  interval  for  recording  data 
extends  for  10  seconds,  from  305.0  to  315.0  seconds.  Within  this  time  frame,  steps  are  made  in 
h  (altitude),  (f>  (roll  angle),  and  /?  (angle  of  sideslip).  The  overlapping,  but  not  simultaneous,  steps 
are  shown  in  Table  A.L  The  lengths  and  timing  of  the  steps  were  set  arbitrarily  with  an  attempt 
to  provide  excitation  in  all  three  axes.  For  identifications  involving  the  base  altitude  of  39,800  feet, 
h  is  set  to  39,850.0  ft,  39,750.0  ft,  and  39,800.0  ft,  respectively  in  the  interval.  Data  are  recorded 
in  the  10  second  interval  for  either  the  unfailed  aircraft  or  the  failed  aircraft  under  any  failure  type 
with  the  failure  occurring  at  305.0  seconds. 

The  identification  techniques  are  the  same  using  the  inverse  identification  formulation  and  the 
indirect  identification  formulation  described  in  Section  3.1.  Differences  appear  in  the  definitions  of 
w[7z]  and  z[n]. 


A. 2  State-Space  Matrices 


For  the  indirect  identification  method,  the  variables  which  comprise  the  w[n]  vector,  given  by 


w[n]  =[ot  q  p  p  T  6hc  ^Ac 

(A.l) 

and  z[n],  given  by 

z[n]  =  [  d  9  ^  p  f 

!  ^ 

i 

(A.2) 

are  recorded  during  the  10  second  interval.  A  least-squares  identification 
matrix  9  is  produced  and  is  given  by 

is  performed  where  a 

6  =  (zw^)(ww^)  ^ 

(A.3) 

where  and  ia  each  w[n]  have  been  filtered  as  described  in  Section  3.2.  The  resulting 

n  X  (n  +  m  +  1)  ^  matrix  contains  the  sub-matrices  A,  B,  and  d. 

Tables  A. 2  and  A. 3  give  the  state-space  matrices  A,  B,  and  d  for  the  unfailed  and  locked 
stabilator  aircraft  with  and  without  differential  stabilator  movement,  respectively,  for  multiple 
flight  conditions.  Flight  characteristics  for  the  unfailed  aircraft  are  given  in  Table  A,4. 
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Table  A.l:  Identification  Stepping  Sequence 


A. 3  Controller  Matrices 

The  controller  matrices  Co,  Go  and  vo  are  obtained  using  (2.15)  and  the  A,  B,  and  d  matrices 
calculated  as  in  Appendix  A. 2,  a-i^d  are  defined  by  (2.10)  and  (2.11),  respectively.  C  is 

defined  as 

‘  0  1  0  0  0  ‘ 

C=  0  0  0  1  0  (A.4) 

0  0  0  0  1  J 

in  order  to  define  y[n]  =  9  P  ^  ]  • 

Tables  A. 5  and  A. 6  give  the  controller  matrices  Cq,  Go,  and  vo  for  the  unfailed  and  locked 
stabilator  aircraft  with  and  without  differential  stabilator  movement,  respectively,  for  multiple 
flight  conditions.  The  matrices  are  obtained  from  the  state-space  matrices  given  in  Tables  A. 2 
and  A. 3. 
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Table  A.2:  State-Space  Matrices  with  doc  — 


a 

9800  ft 

Mach  0.5 

unfailed 

e 

9800  ft 

Mach  0.5 

locked 

b  - 

9800  ft 

Mach  0.9 

unfailed 

f  - 

9800  ft 

Mach  0.9 

locked 

c 

39800  ft 

Mach  0.6 

unfailed 

g  - 

39800  ft 

Mach  0.6 

locked 

d  - 

39800  ft 

Mach  1.4 

unfailed 

h  - 

39800  ft 

Mach  1.4 

locked 
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Table  A. 3;  State-Space  Matrices  with  6oc  —  0 


B 

d 

-0.7599 

1.0026 

0.0244 

|HRir|ua 

-0.0700 

0.0002 

-0.0047 

3.2790 

-6.2749 

-2.9616 

0.1042 

-9.9846 

-0.1242 

0.0179 

-0,0066 

a 

-0.0865 

0.1007 

-0.0902 

mm 

0.1608 

-0.0913 

-0.0611 

0,8685 

-0.1270 

0.4762 

-30.2641 

-2,4635 

mm 

0.8894 

10,7071 

-1.6015 

3.1396 

-0.0331 

0.0230 

4.8450 

0.0473 

0.1660 

-3.0015 

0.1418 

-1.5731 

-0.0030 

-0.0688 

0.0059 

0.0014 

1.9234 

-14.7169 

-3,7149 

0,0538 

-0.0114 

-32.8112 

-0.0703 

-0.0710 

-0.2282 

b 

\  m 

0.0440 

-1.0077 

-0.0631 

-0.0922 

-0.1152 

-0.0735 

-84,4544 

-3.8654 

-0.2023 

0.4927 

22.5359 

-3.6291 

-0.2469 

-0.0217 

20.3607 

-0.0609 

-1.0410 

-0.0938 

0.4239 

-8.9232 

-0.3524 

■1 

-0.2514 

0.9975 

0.0328 

-0.0001 

0.0544 

-0.0399 

0.0089 

0.0111 

2.5802 

■ 

-1.9943 

-0.9338 

-0.2792 

-0,0302 

-0.1843 

-2.8756 

-0.0149 

0.0210 

2.5012 

H 

-0.0985 

0.1124 

-0.1928 

0.2147 

-1.3929 

0.1640 

-0.0245 

-0.0073 

1.9538 

■ 

-0.4102 

0.3859 

-12.4031 

-0.7302 

0.0991 

0.3119 

2.7321 

-0.1330 

6.5433 

■1 

-0,0503 

0.0371 

1.2192 

0.0150 

-0.3652 

0.0196 

0.0140 

-1.1276 

0.3918 

-0.4962 

1,0167 

0,0134 

-0.0045 

0.0044 

0.0190 

-0.0103 

0.0070 

1.2472 

-33.0930 

-1.3887 

0.3751 

0.0251 

0,1680 

-13.3610 

-0.0379 

0.0732 

60.1684 

d 

-0.0008 

0.0347 

0.0596 

0.0501 

-1.0181 

0.0619 

-0.0561 

0.0215 

0.0269 

-0.5169 

-0.0246 

-33.0118 

-1.8412 

0.0058 

-0.1840 

2.8938 

-1.3813 

0.9318 

0.0544 

0.0021 

7.1899 

-0.0474 

-0.5679 

0.0461 

0.0343 

-3.5316 

-0.6290 

-0.7566 

1.0146 

0.0059 

-0.0015 

-0.0098 

-0.0294 

0.0020  . 

-0.0038 

3.3819 

-6.3823 

-2.8453 

0.0109 

-0.0112 

-0.0141 

-4.9151 

-0.0154 

-0.0038 

15.0893 

e 

-0,0449 

-0.1370 

0.3714 

0.1396 

-0.9621 

0.1197 

-0.1602 

-0.0781 

0.5566 

0.1352 

-0.1132 

-29.2927 

-2,4059 

0.7737 

-10.5237 

10.4894 

-1.5874 

-30.8504 

-0.0502 

-0.0388 

4.9104 

-0.0886 

-0.6048 

-1.3384 

0.1439 

-2.9833 

-3.7653 

-1,5705 

1,0129 

0,0051 

-0.0020 

-0.0486 

0,0010 

0.0029 

1.9337 

-14.8281 

-3.7516 

0.4363 

0.0324 

BSRRI 

-16,4431 

-0.0161 

-0.0501 

9.1773 

f 

0.0086 

-0,0545 

-0.1006 

0.0327 

-1.0424 

-0.0042 

-0.0638 

-0.1013 

0.0034 

-0.1417 

-0.0875 

-85,8416 

-3.9473 

-0.2669 

-35.1505 

22.4301 

-3.7710 

-19.5522 

0.0244 

0.2116 

19,9857 

-0.0808 

-3.8319 

0.4734 

-8.9697 

-2.5368 

-0.2671 

0.9858 

0.0387 

0.0010 

0.0366 

-0.0229 

0.0035 

0.0075 

2.8774 

-1.9202 

-0.8826 

-0.1273 

-0.0187 

-0.1196 

-1.4418 

0.0101 

0.0032 

11.7995 

g 

-0.0350 

-0.0600 

0.2628 

0,2281 

-0.8504 

0.1000 

-0.0425 

-0.0127 

0.8804 

-0.1382 

0.0741 

-11.7403 

-0,7709 

0.8784 

-3.5560 

2.5357 

-0.0909 

-23.7294 

-0.0097 

-0.0539 

1.2805 

0.0087 

-0,2809 

-0.5613 

-0.0157 

-1.1260 

-4.1603 

-0.4989 

1.0168 

0.0552 

0.0037 

0.0183 

-0.0066 

0.0023 

1.2489 

-32.5991 

-1.5689 

-0.2683 

-0.0286 

-0.1756 

-6,8991 

0.0337 

-0.0960 

69.0377 

h 

-0.6259 

0.0887 

0.5004 

0.1013 

-0.8338 

0.3101 

-0.0138 

1,9309 

-5.4085 

-0.4432 

-33.0980 

-1.8378 

-0.1025 

-20.0382 

2.9741 

-1.4917 

-0.1756 

-0.1540 

7,4708 

-0.0218 

-0.4275 

-1.6852 

0.0170 

-2,7667 

a 

Mach  0.5 

unfailed 

e  - 

9800  ft 

Mach  0.5 

locked 

b  - 

Mach  0.9 

unfailed 

f  - 

9800  ft 

Mach  0.9 

locked 

c 

Mach  0.6 

unfailed 

g  - 

39800  ft 

Mach  0.6 

locked 

d  - 

Mach  1.4 

unfailed 

ii  - 

39800  ft 

Mach  1.4 

locked 
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Table  A.4:  Flight  Condition  Characteristics 


Altitude 

h 

(ft) 

Speed 

M 

(Mach) 

Airspeed 

(ft/s) 

Air  Density 

P 

(10“^  slug/ft) 

Trim  a 

^trim 

D 

Trim  Sff 
Antrim 

n 

WSSM 

539 

1.8 

4.7 

-2.9 

Bil 

970 

1.8 

1.3 

-0.6 

■B 

581 

0.59 

11.4 

-7.0 

mi 

1355 

0.59 

2.5 

-1.6 

Table  A. 5:  Controller  Matrices  with  Sdc  —  3<5/ic 


Co 

Go 

Vo 

IgilfelsitlB 

-0.6289 

-0.0467 

0.0016 

a 

0.0375 

-0.0244 

1.1633 

-0.0020 

BE  31 

-0.1743 

-0.0012 

-0.0123 

2.5488 

-0.0327 

0.6391 

-0.0433 

-0.4478 

-0.0370 

-0.0028 

0.0013 

-0.0024 

-0.0076 

b 

-0.0052 

-0.0027 

1.1768 

0.0163 

0.0080 

0.0046 

■ 

-0.0073 

-0.0040 

3.2359 

0.0075 

0.1735 

-0.0301 

■ 

0.0025 

0.0051 

-0.6789 

0.5183 

-0.0100 

-0.0070 

-0.0304 

0.7107 

H 

0.2750 

0.0328 

0.0173 

wnujygi 

1.3497 

-0.1835 

-0.1285 

-0.1890 

■ 

0.2138 

-2,2294 

0.0041 

2,2091 

-0.1447 

1.9496 

-0.2820 

II 

-0.0011 

-0.0000 

-2.4785 

0.0012 

0.0002 

4.5061 

y 

0.0744 

-0.0252 

0,0319 

1.0780 

-0.0164 

0.0351 

y 

0.0506 

-0.7605 

0.0184 

3,0158 

-0.0134 

0.6667 

-1.0063 

0.2047 

-0,1086 

0,0229 

2,5436 

-0.0195 

-0.0037 

8.0984 

e 

-0.5871 

0.2504 

-0.1145 

-1.3371 

0.1205 

3.0411 

-0.0272 

-0.0118 

6.0807 

0.1964 

0.0136 

-0.8361 

-0.0031 

1.7783 

-0.0296 

0.6121 

-2.5735 

■ 

-0.3715 

0.0966 

-0,0277 

-2,0737 

-0.1126 

3.3953 

0.0432 

-0.0155 

1.9843 

y 

-0.2669 

0.1205 

-0.0335 

-1.4158 

-0.0278 

4,0984 

0.0463 

-0.0329 

1.7335 

y 

0.0436 

0.0075 

-0,2646 

0.3493 

0.0177 

2.4329 

-0.0000 

0.1443 

-0.4722 

■ 

-4.0665 

0.8917 

-0,1153 

-3.1775 

2,8006 

4.2146 

-0.6431 

-0.2166 

28.6192 

y 

-2.7874 

1,0608 

-0.1185 

-2.1481 

2.0169 

4.9672 

-0,7557 

-0.2429 

23.9967 

y 

0.9281 

-0.0323 

-2.1657 

0.7755 

-0.6959 

1.0446 

0.0459 

1.8954 

-8.9315 

■ 

-2.4374 

-0,2621 

-30.5737 

1.0219 

10.2910 

-0.1971 

0.1248 

72,4688 

u 

-2.4488 

-30,5988 

1.0570 

12.3147 

-0,2309 

0.2084 

74.0099 

y 

SBSSSSS^S 

2.1535 

-0.0792 

2.2263 

0.0038 

0.6619 

-5.3703 

a 

9800  ft 

Mach  0.5 

unfailed 

e 

9800  ft 

Mach  0.5 

locked 

b  - 

9800  ft 

Mach  0.9 

unfailed 

f  - 

9800  ft 

Mach  0.9 

locked 

c 

39800  ft 

Mach  0.6 

unfailed 

g  - 

39800  ft 

Mach  0.6 

locked 

d  - 

39800  ft 

Mach  1,4 

unfailed 

h  - 

39800  ft 

Mach  1.4 

locked 
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Table  A. 6:  Controller  Matrices  with  =  0 


-0.8700 

0.0986 

-0.0139 

-0.1871 

-0.0131 

-0.0026 


-0.0047 

0.9161 

0.0113 

-0.0024 

0.8678 

0.0084 


-0.1520 

-0.2293 

0.0528 

-1.7488 

-2.4201 

0.9055 


-0.3767 

-2.4525 

0.1709 


-0.0001 

0.1123 

0.0060 

0.0069 

0.9950 

-0.0173 


004 

870 

002 


.0001 

.1256 

.8399 


.0007 

.0455 

.2823 


-0.0156 

-0.1062 

-2.2186 

-0.0029 

-0.3397 

-0.7112 


0.0009 

-0.0459 

-0.2815 

-0.0056 

-0.0872 

-2.2163 


.0086 

.3062 

.7261 


.0008 

0.0011 

.2394 

-0.1268 

.0119 

-0.8446 

0.5460 

-0.2016 

0.0399 

0.0831 

0.0146 

0.0018 


86 


-0.1127 

4.6079 

1.1366 

0.0046 

12.4377 

2.1569 


-0.0071 

-0.6467 

0.0051 

0.0025 

-0.2350 

-0.0157 


-0.9014 

-1.3531 

0.3164 

-1.3428 

-1.8041 

0.6859 


0.0762 

0.1070 

0.0505 

1.1311 

1.5343 

0.6331 


0.0149 

4.2618 

2.4467 

-0.0537 

4.5941 

1.0998 


3 


0024 
0460 
6237 
-0.0008 
0.0356 
0.1652 


-0.0506 

0.0617 

1.8930 

0.0148 

0.2614 

0.5498 


0.0019 

0.0665 

-0.0063 

-0.0178 

-0.7060 

0.0265 


-0.0053 

-0.2612 

-0.0050 


-0.0069 

0,0048 

-0.0392 


0.8851 

-2.4798 

0.3319 

4.5029 

-0.0931 

-0,1202 


.0034 

3.0542 

.0193 

5.6481 

.6377 

-2.3599 

-0.0011 

0.0373 

0.1616 

-0.0812 

-0.3880 

2,0167 


-0.0334 

0.1172 

0.6140 


0.5578 

1,6731 

-0.4328 

8.3106 

20.7214 

-8.1265 


10.4436 

73.7459 

-5,5069 


a 

9800  ft 

Mach  0.5 

unfailed 

M 

-  9800  ft 

Mach  0.5 

locked 

b  - 

9800  ft 

Mach  0.9 

unfailed 

R 

-  9800  ft 

Mach  0.9 

locked 

c 

39800  ft 

Mach  0.6 

unfailed 

R 

;  -  39800  ft 

Mach  0.6 

locked 

d  - 

39800  ft 

Mach  1.4 

unfailed 

E 

;  -  39800  ft 

Mach  1.4 

locked 
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Appendix  B 

Calculation  of  the  Local  Angle  of 
Attack  for  Missing  Stabilator 


The  local  angle  of  attack  for  a  missing  stabilator  is  calculated  using  the  following  procedure.  First, 
the  rate  of  change  of  downwash  angle  with  angle  of  attack  is  calculated  by  the  following  equation 
(see  [23]): 

de  2 
da  ~ 

Using  the  parameters  from  Table  2.1  of  5  =  608  ft^  and  h  =  42.8  ft,  tlie  aspect  ratio  is  calculated 
as 


5  608  ft^ 


(B.2) 


Table  B.l  lists  the  coefficients  of  lift  for  various  values  of  angle  of  attack  (a)  and  Mach  number  ( M) 
obtained  from  the  tables  of  the  simulation  code.  The  table  values  are  interpolated  for  the  condition 
of  Mach  0.5.  For  an  8°  change  in  a. 


ACz,  =  0.50  -  (-0.03)  =  0..53  for  Mach  0.4  (B.3) 

ACf,  =  0.52  —  (—0.05)  =  0.57  for  Mach  0.6  (B.4) 

Since  the  flight  condition  of  Mach  0.5  is  desired,  average  the  values  of  ACl  from  (B.3)  and  (B.4) 
to  obtain  ACi  =  0.55.  Calculate 


ACl 

Aa 


0.55  180 

8  TT 


=  3.94=Cl„. 


(B.5) 


Table  B.l:  Coefficients  of  Lift  for  Mach  0.4  and  0.6 


Mach  0.4 

Mach  0.6 

0 

-0.03 

-0.05 

4 

0.24 

0.23 

8 

0.50 

0.52 
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(B.6) 


Substituting  the  values  for  AR  and  in  (B-l)  yields 


da 


2 

7r3.01 


3.94 


0.83. 


The  downwash  angle  e  can  be  calculated  as 


dc 

e  —  Co  + 

da 

The  local  angle  of  attack  at  the  tail  ai  is  then 


(B.7) 


aL  =  a  —  € 

=  a  -  Co  -  0.83a 

=  0.17a -Co  (B.8) 

Unfortunately,  the  constant  cq  cannot  be  determined  easily,  even  using  the  simulation  code.  How¬ 
ever,  it  can  be  viewed  as  a  bias  with  an  effect  similar  to  other  trim  disturbances.  For  simplicity, 
it  was  assumed  that  cq  =  0,  and  the  approximate  local  angle  of  attack  for  the  missing  stabilator 
failure  for  a  flight  condition  of  9,800  feet,  Mach  0.5  was  set  to  0.17a. 
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Appendix  C 

Simulation  Details 


Table  C.l  gives  specific  parameters  for  the  simulations  shown  in  the  figures.  The  failure  using  the 
BLS  identification  occurs  at  305  seconds.  All  RLS  simulations  use  A  =  0.99,  a  =  10.0,  and  failures 
which  occur  at  310  seconds. 
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Table  C.l:  Simulation  Details 


Figure 

Number 

RLS/ 

BLS 

Step 

Size^ 

P 

Value 

Integral 

Action 

Noise 

Present 

Failure 

Type^ 

Limit 

Type^ 

Initial 

Matrices'^ 

1.1 

RLS 

1  k  5 

0 

no 

no 

lock  stab 

none 

(a) 

3.1 

RLS 

1 

0 

no 

no 

lock  stab 

none 

(a) 

3.2 

RLS 

5 

0 

no 

no 

miss  stab 

none 

(a) 

3.3 

RLS 

10 

0,0.1 

no 

no 

miss  stab 

none 

(a) 

4.1 

RLS 

1 

0.1 

no 

no 

lock  stab 

none 

(a) 

4.2 

RLS 

10 

0.1 

no 

no 

lock  stab 

none 

(a) 

4.3 

RLS 

10 

0.1 

no 

no 

lock  stab 

none 

(a) 

4.4 

RLS 

10 

0.1 

no 

no 

lock  stab 

all 

(a) 

4.5 

RLS 

10 

0.1 

no 

no 

lock  stab 

all 

(a) 

5.2 

BLS 

1 

0.1 

yes 

no 

lock  stab 

LS 

(e) 

5.3 

RLS 

1 

0.1 

yes 

no 

lock  stab 

LS 

(a) 

5.4 

RLS 

1 

0.1 

yes 

no 

miss  stab 

LS 

(a) 

5.5 

RLS 

10 

0.1 

yes 

no 

lock  stab 

LS 

(a) 

6.1 

RLS 

1 

0.1 

no 

yes 

lock  stab 

LS 

(a) 

6.2 

RLS 

10 

0.1 

no 

yes 

lock  stab 

LS 

(a) 

7.1 

RLS 

10 

0.1 

no 

no 

i  lock  stab 

LS 

(b) 

7.2 

RLS 

10 

0.1 

no 

no 

lock  stab 

LS 

(c) 

7.3 

RLS 

1 

0.1 

yes 

!  no 

lock  stab 

LS 

(d) 

7.4 

RLS 

30 

0.1 

no 

no 

lock  ail 

all 

(a) 

7.5 

RLS 

50 

0.1 

no 

no 

lock  ail 

all 

(a) 

7.6 

RLS 

10 

0.1 

no 

no 

lock  st/ai 

LS 

(a) 

7.7 

RLS 

30 

0.1 

no 

no 

lock  st/ai 

LS 

(a) 

8.1 

BLS 

10 

0 

no 

no 

lock  stab 

none 

(a) 

8.2 

RLS 

1 

0 

no 

no 

none 

none 

(a) 

8.3 

RLS 

10 

0.1 

no 

no 

lock  stab 

SCI 

(a) 

^  in  ®/s 

^  either  locked  left  stabilator,  missing  left  stabilator,  locked  left  aileron,  or  combined  locked  left  stabilator  and  locked 
left  aileron 

a// represents  all  four  methods  described  in  Chapter  4;  LS  represents  the  least-squares  approximation  to  commanded 
accelerations  method;  SCI  represents  the  scaling  of  control  inputs  method 
*  letter  corresponds  to  the  matrices  listed  in  Table  A. 3 
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Appendix  D 

Gaussian  Random  Number 
Generator 


The  noise  added  to  the  system  is  modeled  as  Gaussian  and  white.  Uniformly  distributed  random 
numbers  between  0  and  +1  are  converted  to  a  Gaussian  distribution  with  a  zero-mean  and  a  unit 
variance  using  the  Box-Muller  method  [21].  A  uniformly  distributed  random  number  generator  is 
written  in  FORTRAN  and  generates  a  random  real  number  between  0  and  -fl  [20].  A  seed,  M, 
must  initially  be  supplied  by  the  user  for  the  first  call  of  the  random  number  generator  function, 
but  subsequent  calls  use  the  seed  from  the  previous  execution  of  the  generator.  Upon  calling  the 
random  number  generator,  M  =  M  ■  65539.  Then,  if  M  is  negative,  M  is  recalculated  as 

M  =  {M +  l)  +  2^ -1,  (D.l) 

where  N  is  the  number  of  bits  defining  the  size  of  words  on  the  particular  machine  that  is  executing 
the  code.  For  this  thesis,  N  =  32.  The  random  number  rand  is  calculated  as 

rand  =  M  •  ^ .  (D.2) 

Two  calls  to  the  rand  routine  generate  two  uniformly  distributed  random  variables  rand  between 
0  and  -f-l  to  be  passed  to  the  Gaussian  conversion  routine.  The  numbers  are  first  transformed  to 
the  range  between  -1  and  -f  1  using  the  following  equations: 

Vi  =  2.0  •  (rand  —  0.5), 

V2  =  2.0 -(rand -0.5),  (D.3) 

r  =  vl  +  vl 

where  rand  is  a  call  to  the  uniformly  distributed  random  number  generator  between  0  and  -fl.  The 
r  variable  is  the  radius  squared  of  the  random  point  on  the  (ui,V2)  plane.  The  r  value  is  tested 
so  that  it  is  always  less  than  1.  Then,  the  region  uniformly  covered  by  (ui,U2)  is  a  circle  and  the 
argument  for  the  square  root  is  positive,  as  well. 

The  two  resulting  random  numbers  Vi  and  t;2  are  remapped  to  two  Gaussian  numbers  using  the 
following  equations: 


gstore  =  vy  •  fac, 
gaus  =  V2  ■  fac. 
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(D.4) 


The  value  in  gstore  is  saved  and  used  during  the  next  call  to  the  gaussian  routine.  The  value  in 
gaus  is  the  output  of  the  routine.  The  parameter  with  noise  is  calculated  as 

param  =  param  +  o'  •  gaus,  (D.5) 

where  param  is  any  of  the  parameters  and  a  is  the  corresponding  standard  deviation  shown  in 
Table  6.1. 

The  same  initial  seed  was  used  for  all  runs  including  noise  shown  in  the  thesis.  However, 
simulations  for  other  seeds  were  also  performed,  without  observing  significant  differences. 
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